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1.  Introduction 

It  has  long  been  known  that  surface  fronts  sometimes  move  at 
higher  velocities  than  the  crossfrontal  component  of  the  actual 
surface  wind  would  support.  In  these  cases,  the  crossfrontal 
component  of  the  geostrophic  wind  has  been  demonstrated  to  provide  a 
somewhat  more  satisfactory  explanation  of  surface  frontal  motion. 
However,  frontal  motion  faster  than  the  geostrophic  crossfrontal 
component  would  support  has  also  been  observed.  Past  research  has 
suggested  that  the  lsal lobar ic  component  of  the  agcostrophic  wind  may 
be  responsible  for  this  phenomenon.  However,  to  date,  the 
isallobaric  wind  has  not  been  shown  to  directly  influence  processes 
which  contribute  to  surface  frontal  motion;  specifically, 
isa 1 lobar lea  1 ly  induced  enhancement  oi  crossfrontal  wind  components 
and  apparent  propagation  through  isallobaric  f rontogenet ica 1 
forcing.  This  paper  will  examine  the  role  of  isallobaric  forcing  on 
surface  frontal  motion  and  f ron togenes is  through  its  influence  upon 
these  two  processes  in  three  cases  studies  over  the  Central  United 
States . 

Studies  of  the  factors  involved  in  surface  frontal  motion  have 
yielded  a  wide  variety  of  results.  For  instance,  Brundidgo  (1965)  in 
his  study  of  cold  fronts,  noted  that  the  frontal  speed  was 
approximately  double  the  normal  wind  component  measured  after  frontal 
passage  at  30  feet  (ft)  and  found  better  agreement  in  most  cases  with 
the  normal  wind  component  measured  at  600  ft.  Brundidgo' s  results 
suggest  that  geostrophic  flow  normal  to  the  frontal  surface  might 
provide  better  agreement  with  the  observed  surface  frontal  motion 
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than  the  actual  surface  winds.  Geostraphic  balance  across  a  front, 
was  used  by  Hoskins  and  Bretherton  <1972)  in  their  model  of 
f rontogenesis,  but  they  noted  that  tangential  accelerations  are  also 
important.  Recent  studies  (Hobbs  and  Pearson,  1982;  Bond  and 
Fleagle,  1985)  have  noted  that  in  some  cases  surface  cold  frontal 
motion  is  accurately  determined  by  the  gravity  current  equation; 
however,  Smith  and  Reeder  (1988)  have  questioned  this  apparent 
agreement  of  observed  surface  frontal  speeds  with  the  calculated 
speeds  of  internal  gravity  waves  and  noted  that  except  where  cold 
fronts  become  orographica 1 ly  trapped,  there  is  little  evidence  to 
suggest  that  frontal  speed  is  controlled  by  gravity  waves  rather  than 
by  processes  which  operate  on  the  frontal  scale  itself,  where  the 
ageostrophic  crossfront  circulation  plays  a  central  role.  One 
component  of  the  ageostrophic  wind  was  considered  by  Brunt  and 
Douglas  in  1928.  In  particular,  they  presented  some  evidence  of 
isallobaric  forcing's  impact  on  surface  frontal  motion,  but  due  to  a 
sparse  observation  network  which  did  not  allow  the  computation  of 
accurate  frontal  positions,  geostrophic  winds,  or  isallobaric  winds 
their  hypothesis  remained  unproved. 

Another  mechanism  by  which  fronts  may  propagate  is  by  the 
f rontogenet lea  1  forcing  associated  with  the  actual  wind. 
Two-dimensional  models  of  horizontal  shear  f rontogenes is  have 
demonstrated  that  when  a  significant  a long- front  temperature  gradient 
exists,  it  is  possible  for  a  cold  front  to  propagate  faster  than  the 
wind  behind  the  cold  front  (Reeder  and  Smith  1986,  1987;  Reeder 

1986).  Furthermore,  Saucier  (1955)  noted  that  the  isallobaric  wind 
Is  f rontogenet lea  1  near  the  regions  of  pressure  falls  and  that  the 
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magnitude  of  the  frontogenesis  in  the  lower  troposphere  is  large  in 
the  same  areas.  In  conditions  similar  t;o  those  noted  above  (strong 
along-front  temperature  gradient),  isnl lobar ic  forcing  on 
frontogenesis  may  also  influence  surface  frontal  motion  through  an 
"apparent’’  propagation. 

In  this  study  the  isallobaric  contribution  to  surface  frontal 
motion  and  surface  frontogenesis  Is  estimated  by  use  of  the 
f r let ion  less  isallobaric  wind.  A  representation  of  typical  pressure 
tendencies  associated  with  a  mid-latitude  cyclone  is  depicted  in  Fig. 

1.  The  frictionless  isallobaric  wind  should  reach  its  maximum 
magnitude  between  the  region  of  pressure  rises  behind  the  cold  front 
and  the  pressure  fall  region  ahead  of  the  cold  front.  Also,  the 
frictionless  isallobaric  wind  should  have  a  local  maximum  of 
magnitude  around  the  pressure  fall  center  north  of  the  warm  front. 
Based  upon  these  distributions,  it  is  expected  that  the  frictionless 
isallobaric  wind  will  be  found  to  contribute  to  increased  speed  of 
surface  frontal  propagation.  Defining  frontogenesis  as  the 
Lagrangian  rate  of  change  of  the  magnitude  of  the  temperature 
gradient,  the  frictionless  isallobaric  wind  in  Fig.  1  would  act  as  a 
frontogenetic  wind  field  through  redistribution  of  the  maximum 
temperature  gradient. 

2.  Data  and  Methods 

2. 1  Data 


Local i ons 


of  the  336  surface  observation  stations  used  to  compile 


T0  -taf 


A  typical  isal lobar ic  distribution  for  an  extratropica 1 
cyclone.  A  <  +  represents  a  pressure  rise  center  and  a  <-) 
represents  a  pressure  fall  center.  Arrows  represent  maxima 
of  the  frictionless  isallobaric  wind  and  broken  lines  am 
isotherms. 
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the  data  base  for  this  study  are  indicated  in  Fig.  2.  Data  input 
from  each  observation  location  included  the  following:  temperature, 
dewpoint,  wind  direction,  wind  speed,  and  altimeter  setting.  The 
above  mentioned  elements  from  hourly  and  corrected  surface  airway 
reports  were  manually  entered  into  computer  input  files.  In  addition 
to  possible  error  from  incorrect  disseminat ion  of  the  observed 
elements,  this  methodology  introduced  the  additional  possibility  of 
error  due  to  incorrect  transfer  of  the  data  to  the  data  sheets  and 
incorrect  typing  into  the  input  fiLes.  The  methods  used  to  screen 
the  data  set  for  errors  are  discussed  in  Appendix  A. 

Brundt  and  Douglas  (1928)  and  Saucier  (1955)  both  noted  that 
correct  representation  of  the  isallobaric  wind  requires  the 
evaluation  of  the  pressure  tendency  field  over  the  shortest  time 
interval  possible.  In  the  absence  of  barometric  traces  for  each 
observation  location  and  given  the  meteorological  data  available,  the 
shortest  time  period  for  which  pressure  tendencies  can  be  computed  is 
1-hour.  Saucier  (1955)  alluded  to  the  use  of  3-hour  pressure 
tendencies  for-  determining  the  pressure  change  due  to  synoptic 
influences,  but  he  noted  that  in  the  middle  latitudes  the  effect  of 
the  semi-diurnal  pressure  wave  is  large  enough  to  make  3-hour  local 
pressure  changes  misleading.  Smaller  scale  influences  such  as 
pressure  changes  due  to  thunderstorms  can  also  influence  the 
tendencies.  In  the  absence  of  a  spectral  analysis  to  identify  the 
different  components  of  the  pressure  tendency,  the  easiest  way  to 
smooth  the  influences  discussed  earlier  is  to  use  the  average  change 


over  a  2-hour  period. 


Locations  of  grid  points  and  surface  stations  used  In  this 
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2.2  Methods 

The  actual  wind  at  any  point  is  defined  as  the  sum  of  the 
geostrophic  and  ageostrophic  winds.  Therefore,  it  is  possible  to 
compute  the  ageostrophic  wind  at  that  point  by  subtracting  the 
geostrophic  wind  from  the  actual  wind.  Actual  winds  at  each 
observation  location  were  broken  into  u  and  v  components  and  fitted 
to  a  one  degree  latitude- longitude  grid  using  t;he  Cressman  <1959) 
objective  analysis  scheme  which  is  described  in  Appendix  B.  The  grid 
is  located  over  the  Central  United  States  and  extends  from  47#N  to  30® 
N  and  from  85°W  to  105®V  <Fig.  2>.  The  total  number  of  grid  points 
is  378.  Gridspacing  is  111.  L9  kilometers  (km)  between  grid  points  in 
the  north-south  direction  and  the  distance  between  grid  points  in  the 
east- west  direction  varies  from  75.83  km  to  96.29  km  at  the  northern 
and  southern  grid  boundaries,  respectively.  Surface  geostrophic  wind 
computations  require  an  accurate  representation  of  the  surface 
horizontal  pressure  gradient  force  which  is  easily  obtained  over 
level  terrain;  however,  the  area  under  consideration  is  sloping 
terrain  which  introduces  the  problem  of  eliminating  the  vertical 
variation  of  pressure  with  height  in  order  to  obtain  the  true 
horizontal  pressure  gradient  force. 

2.2. 1  Horizontal  Pressure  Gradient  Force  Computations 

Sangster  <1960,  1967)  among  others,  has  tinted  the  problems  with 

using  station  pressures  reduced  to  sea  Level  tor  the  representation 
of  the  surface  horizontal  pressure  gradient  force  and  subsequent 


computat  ions  of  the  surface  geostrophlc  wind.  Sa  rigs  ter  <.  1 987)  notes 
that  the  plateau  correction,  which  Is  Implicitly  included  in  the 
reduction  of  station  pressures  to  sea  level  for  station  elevations 
above  305  meters  (m),  has  the  effect  of  producing  a  northerly 
geostrophlc  wind  component  in  the  summer  over  the  Great  Plains,  and  a 
southerly  geostrophlc  wind  component  there  in  winter  with  values  of 
around  10  m  s*  for  the  components  under  extreme  conditions.  Sangster 
<.1960,  1987)  has  also  shown  that  it  is  possible  to  compute  the  slope 

of  a  constant  pressure  surface  in  non- iso baric  coordinates  by  use  of 
the  altimeter  correction  method  developed  by  Bellamy  (1945). 

Hereafter,  all  references  to  Bellamy  will  be  from  his  1945  paper. 

The  altimeter  correction  <D>  is  defined  as  follows: 

H  =  Z  -  Zp,  (1) 

where  Z  is  the  elevation  of  the  station  above  sea  level  and  Zp  is  the 
pressure  altitude  of  the  station  pressure.  Following  Bellamy,  the 
value  of  D  can  bo  used  to  approximate  the  height  of  any  constant 
pressure  level.  For  a  given  station  pressure,  the  pressure  altitude 
is  the  height  of  the  observed  station  pressure  in  a  defined  standard 
atmosphere  using  the  hypsometric  equation 

Zp2  =  -Rdg  ^TplnCPa/PO  +  Zpt,  <l2> 

where  Rd  is  the  dry  gas  constant.,  g  is  the  acceleration  due  to 
gravi t y ,  T/»  '-h  e  mean  standard  atmospheric  temperature  of  the 

layer,  P,  is  i.he  pressure  at  level  1,  Zpj  is  the  pressure  altitude  at 


Level  1 ,  Pj 


i.  -'J> 


the  pressure  at  level  2,  and  Zpa  is  the  pressure 
altitude  at  level  2. 

In  this  study  the  National  Advisory  Committee  tor  Aeronautics 
(NACA)  standard  atmosphere  was  used.  The  NACA  standard  atmosphere 
has  the  following  characteristics  (Haiti ner  and  Martin,  1957): 


a.  Assume  dry  ait  . 

b.  Base  level  pressure,  F'^  -  101:3.25  mb. 

c.  Base  temperature  T  ~  15°C. 

d.  Lapse  rate  -  0.65°C  100  m*  for  z  ^  10759  m. 

e.  gravity  -  980.565  cm  sec  . 


The  correction  for  temperature  and  moisture  may  be  derived  by 
starting  with  the  hypsometric  equation  for  the  actual  atmosphere. 


Z2  -  Zt  =  -Rdg  ‘TvlntPz/Pi > , 


where  Tv  is  the  mean  virtual  temperature  for  the  layer.  Multiplying 

•*  ••  •  / 

Eq.  <3)  by  Tp  Tp  ,  defining  the  pressure  at  level  1  as  1013.215  mb, 
and  substituting  Rq .  <2)  for  the  pressure  altitude  of  the  observed 

station  pressure  <Zp;j)  results  in 


Z2 


ZlO!3. 25  =  TvTp*Zp2, 


where  Zp^  was  calculated  using  Eq .  <5)  Olaltinor  and  Martin,  1957), 


( 

]  . 


ZP2  Cm)  =  44 , 308 1  1 


(  Pz/1  013. 23  ) 


O.  1  902  3 


Now  by  taking  Tv  Tp  as  approx t ma to ] y  equivalent  to  Tv  Tp  wo  are 
assuming  that  the  virtual  temperature  and  the  standard  atmospheric 
temperature  at  the  station  elevation  are  representative  of  the  mean 
temperatures  of  the  column,  respectively.  Following  Haltiner  and 
Martin  (1957),  the  quantity 


Zl013.  25  =  7.7.  ~  TvTp1Zp2> 


represents  the  height  of  the  1013.25  mb  surface.  Bellamy  defined  a 
quantity  called  the  specific  temperature  anomaly  S  which  is  the 
deviation  of  the  actual  virtual  temperature,  corrected  for  gravity 
variations  between  the  actual  and  standard  atmospheres,  from  the 
standard  atmospheric  temperature 


S*  =  (T*-  T p > T p  1 


Bellamy's  definition  of  virtual  temperature  <Tr)  is  given  as 


T*  =  gpg-1Tv, 


where  g f  is  the  acceleration  of  gravity  In  the  standard  atmosphere. 

The  method  used  to  calculate  Tv  and  variations  in  gravity  Is 

w.  -  / 

discussed  in  Appendix  C.  Substitution  of  1  +  S  for-  Tv  Tp  and  D 


for  -  Zpa,  in  Eq .  (5)  yield' 


ZlOl3.  2  f>  =  D  2  S*Zp7. 
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Equation  <9 >  shows  that  the  height  of  the  1013.25  mb  surface  is  the 
sum  of  the  deviation  of  the  observed  column  thickness  from  the 
standard  column  thickness  and  the  correction  due  to  differences  of 
the  observed  virtual  temperature  from  the  standard  atmospheric 
temperature.  Sangster  <1967)  illustrated  that  the  temperature 
averaging  process  used  in  the  reduction  of  station  pressures  to  sea 
level  effectively  "masked"  the  diurnal  variation  of  the  surface 
geostrophic  wind  over  the  Great  Plains.  More  importantly  for  the 
purposes  of  the  present  study  is  the  effect  that  12  hour  temperature 
averaging  would  have  on  an  accurate  determination  of  the  surface 
geostrophic  wind  where  a  change  of  air  mass  has  occurred  due  to 
frontal  passage  since  the  temperature  change  associated  with  some 
fronts  can  be  greater  than  the  normal  diurnal  temperature  change. 

The  slope  of  the  constant  pressure  surface  may  be  obtained  by 
substitution  of  D  in  the  coordinate  transfer  equation 


VnD 


♦ 

VpD  +  (o'Dc'Zp 


1)VnZp, 


(10) 


vhere  the  H  subscript  means  horizontal  in  cartesian  coordinates  and 

-I  M 

the  p  subscript  implies  isobaric  coordinates.  Since  SD^Zp  =  S  ,  the 
hydrostatic  equation  in  the  I)  system,  substitution  of  S  into  Eq. 

<10)  with  some  rearrangement  yields 

*  ♦  *♦ 

7pD  =  VhI)  -  S  VhZp  .  <11) 


Further  examination  of  Eq .  <11)  will  demonstrate  the  problems  of  this 

method  that  Sangster  <1960,  1987)  discusses.  Re i ntroduct ion  of  the 
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subscripts  for  the  pressure  Levels  defined  earlier  and  remembering 
that  the  pressure  altitude  of  the  base  level  (Zpj )  is  equal  to  zero, 

Eq.  (11)  can  be  written  in  the  following  form: 

>  *  *  *♦ 

VpZlOi3.25  =  Cl  +  S  >VhD2  -  S  VhZ2.  (12) 

The  second  term  on  the  right  hand  side  shows  that  the  slope  of  the 
1013.25  mb  surface  depends  partially  on  the  sloped  terrain.  Values 
of  D  computed  at  the  observation  locations  depend  upon  the  station 
elevations  and  if  the  elevations  at  the  grid  points  differ  greatly 
from  the  station  elevations,  introduction  of  a  fictitious  effect 
occurs.  A  satisfactory  solution  to  this  problem  was  achieved  by 
fitting  station  elevations  to  the  grid  using  the  same  Cressman  (1959) 
weighting  scheme  as  was  used  for  computing  the  D  values  and  specific 
temperature  anomalies  at  the  grid  points.  Sharp  changes  in  grid 
point  elevations  can  also  introduce  a  fictitious  slope  that  are  not 
the  result  of  the  horizontal  pressure  differences.  To  minimize  sharp 
elevation  changes,  grid  point  elevations  where  then  smoothed  by  a 
two-pass  five-point  formula  given  in  Appendix  C.  The  smoothed 
elevations  used  are  shown  In  Fig  3.  The  other  problem  results  from 
the  assumption  of  the  mean  standard  atmospheric  temperature  of  the 
column  as  equivalent  to  the  standard  atmospheric  temperature  at  the 
station  elevation.  By  using  the  standard  atmospheric  temperature  at 
station  elevation,  a  temperature  that  will  always  be  colder  than  the 
mean  standard  atmospheric  temperature  is  used.  Referring  to  Eq .  (7) 

the  specific  temperature  anomaly  would  be  greater  than  If  the  mean 
standard  atmospheric  temperature  was  used.  Reduction  of  the  error 


caused  by  this  assumption  was  achieved  by  the  use  of  an  approximate 
mean  standard  atmospheric  temperature  taken  as  the  average  of  the 


13 


surface  and  station  elevation  standard  atmospheric  temperature. 


2.2.2  Geostrophic  Wind  in  the  D  System 


The  expression  for'  the  geostrophic  wind  in  isobarlc  coordinates  is 


♦ 


-gf~*V vz  X  k. 


(13) 


where  f  is  the  coriolis  parameter-  and  k  is  the  vertical  unit  vector. 
Substitution  of  Eq.  (12)  for  the  slope  of  a  constant  pressure  surface 
in  Eq.  (13)  gives 


♦ 

V«  =  -gf^tCl 


S*>7 D  -  S*VZ]  X  k. 


(14) 


The  u  and  v  components  of  the  geostrophic  wind  (u^,  v^ )  obtained  from 
Eq.  (14)  are  as  follows: 


Ug  =  “gf  ‘fCctoory  1  +  -  S* &Zdy  ~  1  ]  (15a) 


Vg  =  -gf  ‘it-dObhc  1  -  S*c1)ch<  l)  +  S*e/Z<?x  1  1  . 


(15b) 


Equations  (15a)  and  (15b)  are  equivalent  to  the  method  used  by  Bonner 
and  Paegle  (1970)  to  confirm  the  diurnal  variation  of  the  surface 
geostrophic  wind  discovered  by  Sangst.er  (1967).  I)  values  and 
specific  temperature  anomalies  were  computed  at  each  station  and  the 
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Cressraan  (1959)  weighting  scheme  was  used  to  fit  the  values  to  the 
grid  points.  Equations  (15a)  and  (15b)  allowed  the  computation  of 
the  u  and  v  components,  respectively,  of  the  surface  geostro  phic 
winds  at  each  grid  point  (excluding  boundary  grid  points)  using 
centered  finite  differencing  to  approximate  the  derivatives. 

2.2.3  The  Ageostrophic  Wind 


The  vector  form  of  the  equation  for  the  ageostrophic  wind  may  be 
written  as  (Saucier,  1955) 


Vag 


=  k  X  (f 


cnrgca. 


where  t  is  time  and  Vn  is  the  actual  horizontal  wind  vector.  Term  1 
in  Eq.  (16)  represents  the  ageostrophic  wind  due  to  a  changing  height 
distribution  with  time  in  isobaric  coordinates  ( isa 1 lohypsic  wind), 
or  in  cartesian  coordinates  is  the  acceleration  due  to  changes  in 
pressure  gradient  and  specific  volume  over  time.  Term  2  represents 
the  ageostrophic  wind  due  to  the  horizontal  advection  of  the 
geostrophic  wind.  Term  3  represents  the  convective  ageostrophic  wind 
due  to  the  vertical  advection  ol  the  geostrophic  wind.  Term  4  is  the 
component  due  to  the  acceleration  of  the  ageostrophic  wind.  Note 
that  the  frictional  term  was  not  included  in  Eq.  (16).  As  shown  by 
Saucier  (1955),  there  is  a  component  of  the  ageostrophic  wind 
(antitriptic  wind)  due  to  frictional  accelerations. 

The  expression  for  the  Isa  1  lobar  in  wind  was  derived  by  Brunt  and 
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Douglas  (1928)  and  can  be  expressed  as  (Haltiner  and  Martin,  1957) 

*  7*  1 

VLS  =  -of  2Vn(dP&l~l )  ,  (1 

* 

where  V,s  is  the  isallobaric  wind  and  *  Ls  the  specific  volume. 

2.2.4  Isallobaric  Wind  in  the  D  System 

Since  the  surface  geostrophic  wind  was  computed  in  isobaric 
coordinates,  consistency  required  an  expression  for  the  frictionless 
isallobaric  wind  in  the  same  coordinate  system.  With  the  assumption 
of  hydrostatic  equilibrium,  the  equivalency  of  the  isallobaric  wind 
in  cartesian  coordinates  to  the  expression  on  the  left  hand  side  of 
Eq.  (18)  can  be  shown. 

-gf  2<1  +  =  — af  2VH(<5P<?t  1  > .  a. 


The  u  and  v  components  of  the  isallobaric  wind  (uJ9  ,  vtt  )  obtained 
from  the  left  hand  side  of  Eq.  (18)  are  as  follows: 


uls  =  S  f  2ui  +  S*  )  (  OC.  1 )  /ctx )  ] , 


(19a) 


*  ~gf  2  r  (1  +■ 


(19b) 


In  order  to  evaluate  Eqs.  (19a)  and  (19b),  the  change  of  I)  was 
computed  from 
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<<5Do*,_1>t  =  CDt  t  -  Dl+1)/2,  (20 > 

where  the  subscript  t  ts  the  time  In  hours.  D  change  values  from  Eq. 
(20)  and  the  specific  temperature  anomalies  at  each  station  were 
computed  and  the  Cressnan  (1959)  weighting  scheme  was  used  to  fit  the 
values  to  the  grid  points.  Centered  finite  differencing  was  used  to 
evaluate  the  derivatives  in  Eqs.  <19a>  and  (19b).  Young  (1973) 
demonstrated  the  deviation  of  the  isallobaric  wind  in  the  boundary 
layer  to  the  left  due  to  the  effects  of  friction  (Fig.  4).  Thus  the 
frictionless  isallobaric  wind  as  an  estimator  of  the  surface 
isallobaric  wind  will  not  passes  the  deviation  to  the  left  as  well  as 
having  a  somewhat  larger  magnitude. 

2.2.5  Locating  Surface  Fronts  and  Computing  Surface  Frontal 
Velocities 

Objectively  defining  a  front  at  the  surface  is  still  at  best  a 
difficult  procedure.  For  example,  defining  a  front  based  upon  the 
magnitude  of  the  gradient  of  temperature  may  work  well  for  anafronts 
(the  prefixes  'ana'  and  'kata'  are  as  defined  by  Bergeron  (1937)) 
which  are  usually  defined  by  sharp  changes  in  temperatures  however, 
the  katafront  is  usually  accompanied  only  by  a  drop  in  relative 
humidity  (Browning  and  Monk,  1982).  An  example  of  the  magnitude  of 
the  temperature  gradient  for  one  hour  from  the  31  Jan  89  case  study 
illustrates  the  difficulty  of  defining  the  front  based  soley  on 
magnitude  of  the  temperature  gradient  (Fig.  5).  It  is  generally 
accepted  that  the  maximum  temperature  gradient  should  be  located 


Fig.  4. 


Wind  vector  hodographn  for  Eknwn  (Vo)  flow  and  frictional 
isal lobar ic  flow  (Vr  )  .  After  Young,  1973 
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behind  a  cold  front.  If  this  criterion  is  used  as  the  only  parameter 
by  which  to  locate  the  front,  the  western  one-third  of  the  cold  front 
is  well  defined;  however,  the  eastern  two-thirds  of  the  cold  front 
does  not  conform  to  this  pattern.  Based  upon  the  uncertainty  of 
frontal  postioning  in  areas  of  weak  temperature  gradient,  it  was 
decided  to  use  the  more  subjective  approach  of  defining  the  surface 
frontal  positions  using  all  available  surface  parameters.  In  order 
to  reduce  the  bias  that  may  result  from  this  procedure,  surface 
analyses  produced  at  the  National  Meteorological  Center  were  used  as 
guides  for  location  of  the  surface  fronts  at  1500,  1800,  2100,  and 

0000  G. M. T. .  Surface  frontal  positions  at  intermediate  hours  were 
determined  from  plotted  station  data. 

Frontal  velocities  were  determined  from  hourly  surface  positions 
for  various  hours  within  each  case  study.  Some  error  in  the 
determination  of  surface  frontal  velocities  is  to  be  expected  with 
the  subjective  method  of  frontal  location;  however,  the  resultant 
frontal  velocities  are  considered  accurate  enough  for  the  purposes  of 
this  study.  The  component  of  the  geostrophic,  actual,  and 
frictionless  isallobaric  winds  normal  to  the  surface  front  were 
determined  for  the  same  periods  as  the  frontal  velocities.  The 
normal  component  of  the  respective  wind  for  each  segment  of  the 
frontal  zone  (segment  defined  as  the  length  of  front  one-half  the 
distance  from  each  grid  point)  is  determined  from  the  nearest  grid 
point.  The  nearest  grid  point  definition  is  based  upon  one-third  of 
the  distance  from  the  middle  of  the  front  to  the  grid  point.  If  the 
front  was  located  more  than  this  distance  away  from  the  two  nearest 
grid  points,  the  crossfrontal  component  was  taken  as  a  linear  average 
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Fig. 


Fig,  6 


Magnitude  of  the  temperature  grad  ten  t  \\VTU  f  rora  31  January  1989 
far  1900  G.M.T.  Units  are  ®K  Km*  X  108  with  a  contour  of 


Hourly  frontal  positions  over  the  central  United  States  for 
1800-0100  G.M.T.  31  January-1  February  1989. 
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of  the  crossfrontal  components  from  the  two  nearest  grid  points. 
Frontal  velocities  along  with  the  crossfrontal  component  by  the 
actual,  geostrophic,  and  frictionless  isallobaric  winds  are  displayed 
together  in  figures  hereafter  called  frontal  velocity  maps. 

2.6  Computation  of  Surface  Frontogenet ica 1  Forcing 

Petterson  (1056)  demonstrated  that  by  assuming  adiabatic  motion 
and  neglecting  the  vertical  component  of  motion,  the  kinematic 
factors  involved  in  surface  f rontogenesis  of  a  potential  temperature 
field  may  be  expressed  in  terms  of  the  deformation,  divergence,  and 
the  angle  of  the  isotherms  from  the  axis  of  dilitation.  Petterson 
further  demonstrated  that  for  synoptic  scale  motions  the  divergence 
term  may  be  neglected  and  obtained  the  now  classical  result  that  if 
the  angle  from  the  isotherms  to  the  axis  of  dilitation  is  less  than 
45  degrees  f rontogenes is  occurs  and  if  the  angle  exceeds  45  degrees 
frontolysis  occurs.  In  the  scale  of  motion  associated  with  a 
developing  cyclone  or  a  front,  it  is  possible  for  the  divergence  term 
of  the  surface  f rontagene t ica 1  function  to  be  of  the  same  order  of 
magnitude  or  larger  than  the  deformation  term.  For  example,  in  their 
case  study  of  a  developing  cyclone,  Moore  and  Blakey  <1988)  had 
contributions  to  the  surface  f rontogenesis  from  the  divergence  term 
of  the  same  order  of  magnitude  as  the  deformation  term.  By 
definition  the  contribution  to  surface  f rontogenesis  by  convergent  or 
divergent  motions  of  the  actual  wind  field  has  to  come  from  the 
ageostrophic  component.  Surface  i rontogenes i s  evaluated  as  the 
Lagrangian  rate  of  change  of  the  horizontal  temperature  gradient 
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magnitude  and  with  the  vertical  component  of  motion  removed  may  be 
written  as  (Jascourt  et  al. ,  1988): 

$  =  d(  IIVTII)/db  =  -i cfYctK ~ 1 1  ‘ HAicfK  1+  ffldy  ll 

+  cfUty  1  i  <7T<?x  1  <9u<?y  1  +  cO'c^y  "  1  <?v<?y  "  1 1  - 
[  dTcb<  4<?CdTdt  1  >/?x  1  +  cfldy  “  ‘d'CdTdt  1  ><?y  11>/IIVTII, 


where  llVTll  is  the  magnitude  of  the  temperature  gradient.  By  using 
the  quantities  defined  below 


Shearing  Deformation  =  Sr  =  dfudty  1  +  (Jvcfx  ~ 1 , 

Stretching  Deformation  =  St  =  OurJx ~  1  -  OvOy  ~ 1 , 

♦  ♦ 

Divergence  =  V  •  V  =  Div  =  <?uc>x  +  <?v<?y  1 , 


it  can  be  shown  that  Eq .  <21)  may  be  written  as 

8  =  -iCcrrcK^cn'dy  1  >sr  +  [(<da<?x1)2+(cyi’<yy”1)2>Div]/2 

+  CtCcO'cbc-1  )2-fca’c)y“1>2)St  1/2  -  <22> 

[dTdbc  "1^ CdTdt  _l)dtx  1  +  c^Toy  ^CdTdt " 1  ><?y ~ 1  ]  >/ll 7TII . 

The  last  term  on  the  right  hand  side  of  Eq.  (.22)  represents  the 
temperature  tendencies  portion  of  the  f rontogenesis  function 

icfl&K  ‘^C.dTdt-  1  )<?x~ 1  +  dTcJy"  1<XdTdt  "  1  )c^y  ~ 1  J/llWlI .  <23) 

Other  studies  (i.e. ,  Jascourt  et.  al. ,  1988)  have  often  ignored 

this  term.  A  short  discussion  on  the  validity  of  neglecting  the 


temperature  tendencies  term  Is  presented  in  the  first  case  study. 


Breaking  Eq.  (23)  into  its  advectlve  and  local  change  components  and 
reuniting  the  result  with  Eq.  (22)  gives  a  form  of  the 
f rontogenet ica 1  function  that  can  be  evaluated  far  the  u  and  v 
components  of  the  actual  and  isal lobar ic  winds, 

g  s  -{[Cd'T<?x_1<9T<?y  4>Sr]  +  l<(<3T<?x  1 ) 2  +  (cfYOy  1>z)Div]/2 
+  [  C  C<9T<9x  1 )  2  -  (cTl'cty  1)2>St  1/2  -  L  cTTcfx  "‘cKOTdt ' 1  >&x  ~  1 

+  <5rT<?x_  4<X  V»VT)d>x_  1  +  crt’dty  ~  ~ 1  ><?y  1  <24 

+  cfldy  ~  V»VT)<?y  1  >211 VI’  II , 

♦  * 

where  V«VT  was  evaluated  as 

V«Vr  =  uCdTVcbO  +  vCcOVtfy)  . 

Equation  <24)  was  evaluated  for  the  actual  and  frictionless 
isallobaric  winds.  The  frictionless  isallobaric  wind  was  available 
on  all  except  the  exterior  grid  points  (latitudes  47dN  and  30*N, 
longitudes  105®W  and  85® V).  The  derivatives  in  Eq.  (24)  were 
approximated  using  centered  finite  differencing  except  for  the 
exterior  grid  points  where  uncentered  finite  differencing  was  used. 
Because  the  times  when  the  actual  wind  is  measured  at  observation 
locations  can  vary  considerably,  a  sensitivity  study  was  performed  to 
determine  how  f rontogenet ica 1  forcing  can  be  influenced  by  small 
changes  in  the  actual  wind,  This  is  discussed  in  Appendix  D. 

3.  Results 


Contri  biit  ions  by  the  isallobaric  wind  to  surface  frontal  motion 
and  surface  f rontognet ica 1  forcing  are  examined  in  three  case 
studies.  Differing  synoptic  conditions  in  each  case  study  allow 


observation  of  varied  isallobaric  patterns  in  the  vicinity  of 
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surface  fronts.  The  31  January  1989  case  study  Investigates  the 
isallobaric  effects  on  a  cold  front  that  serves  as  a  boundary  for  an 
arctic  air  mass  moving  south  from  Canada.  The  other  two  case  studies 
( 15  November  1988  and  7  January  1989 >  examine  isallobaric  influences 
associated  with  developing  mid- lat itude  cyclones. 

3. 1  31  January  1989  Case  Study 

On  31  January  1989  a  cold  front  located  over  South  Dakota  and 
Minnesota  at  i800  G.M.T.  moved  south  and  was  located  from 
south-central  Nebraska  through  Iowa  and  into  central  Wisconsin  by 
0100  G.M.T.  1  February  1989  <Flg.  6).  Referring  to  Fig.  6,  it  is 
apparent  that  the  frontal  zone  did  not  d-  nstrate  a  constant 
movement;  decreased  frontal  sr°e.’  Is  observed  during  the  period 
2000-2100  G.M.T.  fallowed  by  increased  frontal  speed  from  2200-0000 
G.M.T.  These  frontal  velocity  changes  correspond  to  pronounced 
changes  in  the  D  value  (pressure)  tendency  field  during  the  same 
period.  The  National  Meteorological  Center  (NMC)  1200  G.M.T.  850  mb, 
700  mb,  and  500  mb  maps  (Figs.  7,  8,  and  9)  reveal  that  the  southerly 

push  of  cold  air  through  the  period  of  interest  was  in  response  to 
northerly  flow  in  the  lower  levels  created  by  an  intense  closed  law 
aver  Hudsons  Bay  and  the  high  pressure  ridge  over  Northwestern  North 
America.  Closer  scrutiny  of  the  upper  air  maps  reveals  several  other 
features.  For  example,  the  warm  air  over  Nebraska  and  Kansas  is  the 
result  of  flow  over  the  Rocky  Mountains  as  indicated  by  the  lee  side 
trough  in  the  lower  troposphere .  Surface  temperatures  in  response  to 
this  adiabatic  warming  reached  maxi  mums  In  the  mid  to  upper  70' s  over 


except  at  700  mb. 
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Kansas  and  southern  Nebraska.  Notice  the  short  northwesterly  tilt  of 
the  trough  up  to  700  mb  above  which  It  became  almost  vertically 
stacked  to  500  mb  just  north  of  the  United  Sta tes-Canad Ian  Border. 

In  addition,  a  decreasing  temperature  gradient  with  height  is 
observed . 

The  gridded  surface  temperature  field  at  the  beginning  of  the 
case  study  is  shown  in  Fig.  10.  This  remarkable  gradient  created 
some  very  interesting  temperature  and  pressure  changes  as  the  air 
mass  moved  southward.  For  example,  at  2200  G.M.T.  Valentine, 

Nebraska  reported  a  temperature  of  69°F  while  at  Pierre,  South  Dakota 
(located  approximately  150  km  to  the  north)  the  temperature  was  26® 

F.  At  2300  G.M.T. ,  one  hour  after  frontal  passage  occurred, 

Valentine  reported  a  temperature  of  32®F.  During  the  same  one  hour 
period  when  the  temperature  drop  occurred,  Valentine  reported  an 
incredible  pressure  rise  of  6.4  mb!  Precipitation  associated  with 
the  front  was  in  the  form  of  light  snow  well  behind  the  surface 
front.  Intense  pressure  rises  behind  the  front  along  with  pressure 
falls  ahead  of  the  front,  coupled  with  a  frontal  zone  well  defined  by 
temperature,  made  for  an  ideal  case  study  to  test  the  hypothesis. 

The  1900  G.M.T.  actual,  geostrophic,  and  isallobaric  winds  In  the 
vicinity  of  the  frontal  zone  are  indicated  in  Figs.  11-13, 
respectively.  Frontal  segment  labels  and  the  ageostrophic  winds  at 
the  same  grid  points  used  to  compute  the  crossfrontal  components  by 
the  actual,  geost roph ic ,  and  isallobaric  winds  are  displayed  in  Fig. 
14.  The  most  apparent  difference,  illustrated  by  the  1900  G.M.  T. 
winds,  is  the  larger  degree  of  perpendicularity  to  the  frontal  zone 
by  the  actual  and  isallobaric  winds  as  compared  to  the  geostrophic 


X/ 
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Fig  13. 


Same  as  Fig.  11,  bul.  far  isal  lobarla  winds. 
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Fig.  14.  1900  G.R.T.  ageostrophic  winds  plotted  at  the  grid  points 

used  to  determine  the  crossfrontal  components  by  the 
actual,  geostrophic,  and  isallobaric  winds.  Letters 
designate  frontal  segments  on  the  superimposed  front. 


Fig.  15.  Frontal  velocity  map  for  the  period  1900-2000  G.M.T.  31 

January  1989,  Frontal  velocities  Cm  s* )  are  indicated  by 
the  larger  numbers  and  correspond  to  the  frontal  segments 
Indicated  in  parenthesis  by  the  frontal  velocity.  Arrows 
Indicate  the  approximate  direction  of  motion.  Small  numbers 
indicate  the  crossfrontal  component  (m  s"* )  by  the 
aetua  1  /geost.roph ic/ isa  1  lobar ie  winds.  Positive  (negative) 
cross  frontal  components  Indicate  flow  In  the  same 
(opposite)  direction  of  frontal  motion. 
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winds,  especially  over  fronts  I  segments  0-1.  However,  the  larger 
magnitude  of  the  geostrophic  winds  over  frontal  segments  B-D  and  F 
made  for  positive  (in  the  direction  of  motion)  crossfrontal 
components  (Fig.  15)  that  were  generally  similar  to  the  observed 
frontal  motion.  The  1900  G.R.T.  frontal  velocities  map  (Fig.  15) 

shows  that  the  geostrophic  crossfrontal  components  do  not  support  the 
observed  frontal  motion  at  segments  A  and  G- I .  The  actual  wind 
crossfrontal  components  are  much  closer  to  the  observed  velocity  at 
these  segments.  The  reason  may  be  due  to  the  ageostrophic 
crossfrontal  component.  The  question  then  becomes  how  much  of  the 
ageostrophic  crossfrontal  component  is  due  to  the  lsallabaric 
acce lera t ion. 

Assuming  that  the  frictional  and  isallobaric  forces  are  the 
dominant  components  of  the  observed  ageostrophic  wind,  it  is  possible 
to  estimate  the  isallobaric  contribution  by  comparison  of  the 
observed  ageostrophic  wind  to  the  ageostrophic  wind  resulting  from 
friction.  Frictional  effects  wi 1 L  be  parameter i zed  by  using  a  40 
degree  cross-contour-  flow  and  30%  reduction  (a  conservative  value)  of 
the  geostrophic  wind  speed  at  three  grid  points  used  to  compute 
crossfrontal  components.  Starting  with  segment  A  (Fig.  16(1)),  the 
ageostrophic  wind  due  to  the  frictional  acceleration  (Vagf)  has  a 
much  sraalLer  magnitude  and  would  riot  posses  the  normal  component  that 
the  observed  ageostrophic  wind  <Vag)  posses.  The  ageostrophic  wind 
due  to  the  frictional  parameterization  at  segment  D  (Fig.  16(2))  is 
similar  to  the  observed  ageostrophic  wind.  However,  since  the  30% 
reduction  in  goost.roph lc  wind  speed  is  probably  too  small,  the 
ageostr oph 1 c  wind  due  to  friction  would  have  a  larger  magnitude  and 
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Fig.  16.  <a>.  Ageostrophic  wind  <Vagf)  at  1900  6.M.T.  from  estimate 

nf  the  actual  wind  (Vj,  >  dr^e  to  frictional  parameterization 
on  the  geostrophic  wind  <Vg>.  (b>  .  Comparison  of 

the  ageostrophic  wind  due  to  fricition  (Vagf)  with  the 
observed  ageostrophic  wind  (Vag)  and  the  f r ici t ionless 
isal lobaric  wind  (Vis).  Geostrophic  and  ageostrophic  winds 
are  from  the  grid  points  used  to  compute  the  crossfrontal 
components  at  frontal  segments  A  (1),  D  (3),  and  G  (3). 
Vector  magnitude  is  based  upon  a  ratio  of  two  millimeters 
to  one  m  s"1 . 
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point  more  towards  the  northeast.  This  suggests  that  the  isallobaric 
acceleration  is  counteracting  the  frictional  acceleration  resulting 
In  a  larger  actual  wind  crossfrontal  component.  This  simple  analysis 
demonstrates  that  it  is  not  necessary  for  the  ageostrophic  wind  to 
posses  a  positive  crossfrontal  component  to  Indicate  an  isallobaric 
contribution.  Segment  G  (Fig.  16(3))  represents  a  case  where  the 
geostrophic  crossfrontal  component  is  negative,  but  the  positive 
actual  wind  crossfrontal  component  is  primarily  the  result  of 
friction.  The  above  examination  suggests  that  the  magnitude  of 
Isallobaric  influence  on  the  actual  wind  crossfrontal  components  can 
vary  from  dominant  to  insignificant. 

Frontogenet i cal  forcing  by  the  actual  wind  and  contributions  by 
the  various  terms  (shearing  and  stretching  deformation,  divergence, 
and  temperature  tendencies)  in  Eq.  (24)  are  shown  in  Figs.  17-21. 
Given  the  number  of  points  that  comprise  the  area  of  interest,  a 
simple  analysis  to  estimate  the  f rontogenet ica 1  forcing  due  to  the 
ageostrophic  wind  resulting  from  friction  is  beyond  the  scope  of  this 
study.  Magnitudes  of  the  component  terms  were  not  smoothed  prior  to 
analysis.  Notice  that  in  areas  where  pronounced  cyclonic  shearing  is 
indicated  (western  South  Dakota  and  northwestern  Jowa),  the 
f rontogenet ica 1  forcing  due  to  shearing  deformation  of  the  actual 
wind  (Fig.  17)  is  contributing  to  frontoiysis.  As  in  Moore  and 
Blakey's  (1988)  results,  the  magnitude  of  the  divergence  term  (Fig. 

19)  is  equivalent  to  the  magnitude  by  the  deformation  terms  (Figs.  17 
and  18).  'the  f  ran  togonos  i  s  maximum  due  to  the  divergence  terra  (Fig. 
19)  along  frontal  section  A-D  occurs  in  an  area  where  the  actual 
winds  (Fig.  11)  are  convergent,  within  with  a  north-south  temperature 
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Fig.  17.  Frontogenetlcal  forcing  due  to  shearing  deformation  of  uhe 
actual  surface  wind  at  1900  G.M.T.  31  January  1989 
(X  10  ®K  km1  hr).  Positive  (.negative)  values  represent 
f rontogenesls  < f rontolysls) .  Surface  frontal  position  is 
super  1 mposed . 
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gradient  (Fig.  LO).  One  of  the  more  Interesting  results  of  this 
study  was  the  significance  of  the  temperature  tendencies  term  (Fig. 

20)  contribution  to  the  total  f rontogenet lea  1  iorcing  by  the  actual 
wind  (Fig.  21).  Notice  that  frontal  segments  G  and  I  (Fig.  21)  are 
in  a  f rontogenet ica 1  farcing  pattern  indicative  of  movement  in 
response  to  a  redistribution  of  the  temperature  gradient  ("apparent" 
propagation) . 

Contributions  by  the  isallobaric  f rontogenet ica 1  forcing  terms 
(Figs.  22-25)  have  somewhat  larger  magnitudes  than  for  the  actual 
wind;  this  is  not  surprising  given  the  larger  magnitudes  of  the 
frict ion  less  isallobaric  wind.  The  convergent/d i vergent  nature  of 
the  frictionless  isallobaric  wind  is  illustrated  by  the  divergence 
term  (Fig.  24).  As  expected,  the  magnitudes  of  the  divergence  term 
are  larger  than  for  either  shearing  (Fig.  22)  or  stretching 
deformation  (Fig.  23).  The  divergence  term  contribution  is  easily 
visualized  by  observing  the  D  value  tendencies  (Fig.  27)  together 
with  the  temperature  gradient  <Fig.  10).  In  addition,  frontal 
section  A-D  and  to  a  lesser  extent  frontal  segments  H  and  I  are 
located  in  a  f rontogenet ica 1  farcing  pattern  due  to  the  divergence 
term  conducive  of  "apparent"  propagat ion .  An  approximate  estimate  of 
the  pressure  change  in  millibars  may  be  obtained  from  Fig.  27  by 
using  the  ratio  of  one  millibar  pressure  change  to  10  meters  of  D 
value  change  as  required  for  hydrostatic  equilibrium  near  the 
surface.  The  temperature  tendencies  term  (Fig.  25)  was  the  primary 
contributor  to  the  total  f rontogone t lea  1  forcing  by  to  the 
Isallobaric  wind  (Fig,  26).  Closer  examination  of  tire  total 
f rontogenet Ica 1  forcing  (Fig.  26)  reveals  a  pattern,  through  which 


Frontogenet ica 1  forcing  due  to  shearing  deformation  of  the 
lsal lobaric  wind  at  J900  G.M.T.  31  January  1989 
(X  Hr  K  km  hr).  Positive  (negative)  values  represent 
f rontogenes is  ( f rontolysis) .  Surface  frontal  position  is 
super i mposed . 


Fig.  26, 


Same  as  Fig.  22,  but  for  tota L  isallobaric  f rontogenet i ca 1 
f  arc i ng . 
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Fig.  27.  Tendency  of  D  values  at.  1900  G.M.T.  31  January  1989  (mhr*  ). 

Height  falls  are  indicated  by  dashed  lines  and  height  rises 
by  solid  lines.  Surface  frontal  position  is  superimposed. 
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frontal  section  F~  l  lies,  favorable  for  "apparent”  propagation.  The 
coincidental  agreement  of  the  isallobaric  and  actual  total 
f rontogenet lea  1  farcing  over  frontal  segments  G  and  1  would  suggest 
that  some  of  the  observed  frontal  motion  was  indeed  occurring  through 
f rontogenet lea  1  forcing  produced  by  the  Isallobaric  wind. 

The  magnitude  of  the  temperature  tendencies  term  has  implications 
for  future  f rontogenet lea  1  forcing  studies.  Argument  for  inclusion 
of  the  term  can  best  be  made  by  examination  of  the  assumption  under 
which  it  has  frequently  been  eliminated.  For  example,  Jascourt  et. 
al.  <1988)  neglects  the  term  based  upon  an  observed  pattern  of  little 
temperature  change  over  time.  Recall  from  Eq.  (24)  that  the 
temperature  tendencies  term  is  a  combination  of  the  gradient  of 
temperature  change  and  the  gradient  of  advection.  Figure  28 
identifies  the  location  of  two  grid  points  ahead  and  behind  the  front 
were  the  contribution  by  the  gradient  of  temperature  change  and  the 
gradient  of  advection  for  each  wind  have  been  calculated  and 
presented  in  Table  1.  It  is  apparent  from  Table  1  that  the 
contribution  to  the  temperatures  tendencies  terra  by  the  gradient  of 
advection  can  be  of  equal  or  greater  magnitude  than  the  contribution 
by  the  gradient  of  local  temperature  change.  These  calculations 
demonstrate  that  even  with  strong  temperature  change,  the  spatial 
gradient  of  advection  can  be  the  dominant  portion  of  the  temperature 
tendencies  term. 

Isallobaric  forcing  has  been  subjectively  demonstrated  as  a 
contributor  to  the  actual  wind  crossfrontal  components  at  frontal 
segments  (A,  D,  H,  and  I).  At  other  frontal  segments  <B,  C,  E,  F, 
and  G>  isallobaric  influence  on  the  actual  wind  crossfrontal 
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components  was  small  or  unclear.  Differences  between  the  observed 
frontal  motion  and  the  actual  wind  crossf ronta 1  components  over 
frontal  segments  G  and  I  are  not  explained  by  the  f rontogenet ica 1 
forcing  by  the  actual  wind.  Magnitudes  of  the  actual  frontogenetical 
forcing  give  a  temperature  gradient  realignment  over  100  km  of  only 
. 5°C  for  segment  I  Increasing  to  1.5°C  for  segment  G .  Given  the 
magnitude  of  the  temperature  gradient  associated  with  the  front  (Fig. 
5),  the  amount  of  "apparent”  propagation  induced  by  isallobarie 
farcing  is  does  not  appear  large  enough  to  account  for  the  difference. 
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Table  1.  Frontogenetical  forcing  by  the  temperature  tendencies  terra 
(Total),  contribution  by  the  gradient  of  temperature  change,  and 
contribution  by  the  gradient  of  advection.  <.X  10^°K  kn/hr1). 

The  period  2200  to  2300  G.M.T.  marked  the  beginning  of  Increased 
frontal  speed  over  frontal  section  A-K  (Figs.  29  and  30).  Note  the 
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strength  of  the  isallobarie  crossf rontal  components  in  response  to  a 
north-south  D  change  gradient  over'  frontal  section  E-G  (Fig.  31). 


This  intense  gradient  was  the  result  of  the  D  rise  center  having  a 
closer  proximity  to  the  front  than  in  previous  hours.  The 
isallabar ic  influence  on  the  actual  wind  crossfrontal  components  over 
frontal  segments  E  and  F  is  unclear-  due  to  friction;  however,  a 
significant  isallobaric  contribution  is  suggested  on  frontal  segment 
G  by  virtue  of  a  positive  ageostrophio  crossfrontal  component  (Fig. 
29)  not  indicative  of  friction.  .Smaller  isallobaric  contributions 
are  possible  over  frontal  segments  A,  G,  H,  I,  and  Q. 

Frontogenet.  ica  1  farcing  due  to  the  actual  and  isallobaric  winds 
(Figs.  32  and  33,  respectively)  during  tills  period  also  demonstrated 
an  interesting  change.  This  change  was  the  f  rontolysis  by  both  winds 
behind  frontal  section  D-F.  With  the  reposition  of  the  D  value  rise 
center  over  the  maximum  temperature  gradient,  the  divergence  term 
became  a  substantial  contributor  to  the  isallobaric  frontolysis. 
Disappointingly,  the  magnitude  of  the  temperature  gradient 
realignment  due  to  the  actual  wind  frontogenet  lea  1  forcing  (Fig.  32) 
is  only  2®C  over  100  km  and  is  not  large  enough  to  explain  the 
difference  between  the  actual  crossfrontal  components  and  the 
observed  frontal  motion.  The  coincidental  occurrence  of  increased 
frontal  velocity  with  on  increased  D  value  rise  center  in  close 
proximity  to  the  front  suggest  an  isallobaric  contribution.  The 
diagnosed  isallobaric  contribution  has  been  subjectively  illustrated 
on  both  the  actual  wind  crossfrontal  components  and  on 
f rontogenet i ca I  forcing.  The  limited  geographical  extent  of 
Influence  relegates  the  isallobaric  role  to  one  of  many  potential 
cont r i tutors. 

By  2300  G.M.T.  a  west. -southwest  shift,  of  a  further  increased  D 
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Location  of  the  points  listed  In  Table  1.  Surface 
frontal  position  Is  superimposed. 
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Fig.  30  Same  as  Fig.  15,  except  for  2300  G.M.T.  . 


Fig.  31. 


Same  as  Fig.  27,  except  for  2200  G.M.T.  . 


Fig.  33.  Same  as  Fig.  26,  except  at  2200  G.M.T. 
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value  rise  center  to  extreme  southern  South  Dakota  (.Fig.  34) 
correspondingly  moved  the  largest  isal lobaric  crossfrontal  components 
through  frontal  section  A-.J  (Figs.  35  and  36).  Notice  that  frontal 
acceleration  (from  the  previous  hour)  is  indicated  over  frontal 
sections  A-F  and  K-L  (Fig.  36).  Counteraction  of  friction  over 
frontal  segments  A,  C,  D,  and  F  is  indicated  by  ageostrophic  win  is 
parallel  to  the  front  (Fig.  35)  and  actual  crossfrontal  components 
similar  to  geostrophic  crossfrontal  components  (Fig.  36).  Unlike 
previous  hours,  a  somewhat  consistent  isal lobaric  influence  was 
indicated  where  the  maximum  isal lobaric  crossfrontal  components 
occurred.  Unfortunately,  no  improvement  of  frontal  motion 
diagnostics  by  the  actual  crossfrontal  components  is  indicated  at  any 
frontal  segment  where  an  lsallabaric  influence  is  suggested. 
Isallobaric  influence  (Fig.  37)  on  the  actual  wind  f rontogenet ical 
forcing  (Fig.  38)  would  again  appear  to  be  manifested  as  frontolysis 
behind  frontal  section  B-L.  While  the  addition  of  velocity  induced 
by  the  actual  wind  f rontogenet ical  forcing  to  the  actual  wind 
crossfrontal  components  would  appear  to  provide  better  agreement  with 
observed  frontal  motion  than  the  geostrophic  crossfrontal  components 
at  some  segments  (C,  D,  and  F) ,  it  is  unclear  that  this  would  be  the 
case  at  other  segments  (A,  B,  and  E).  The  limited  area  of 
isallobaric  influence  is  shown  by  little  or  contribution  on  the  other 
frontal  section  (K-L)  displaying  increased  velocity. 

In  summary,  significant  isallobaric  kinematic  and  f rontogenet ica 1 
forcing  contributions  were  indicated  at  some  frontal  segments  in  this 
case  study.  The  appearance  of  "apparent"  propagation  in  this  case 
study  suggests  that  some  frontal  motion  occurs  through  frontal 
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Fig. 


Fig. 


37.  Same  as 


Fig.  33,  except  at  2300  G.M.T. 
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rede  ve  Lopraen  t .  The  isallobaric  acceleration  has  been  shown  as  a 
contributor  to  this  f rontogene t lea  1 1 y  Induced  motion.  Diagnosing  an 
isallobaric  contribution  t;o  surface  frontal  motion  thi  ough  the  actual 
wind  crossfrontal  components  makes  the  significance  of  isallobaric 
crossfrontal  components  difficult  to  determine.  However,  strong 
isallobaric  crossfrontal  components  were  indicated  across  many 
frontal  segments.  The  fact  that  the  actual  wind  crossfrontal 
components  on  these  segments  did  not  consistently  (until  the 
isallobaric  acceleration  became  unusually  large)  display  an 
isallobaric  influence  suggests  that  the  contribution  may  have  been 
somewhat  random.  When  some  consistency  of  isallobaric  influence  on 
the  actual  wind  crossfrontal  components  was  indicated  (2300  G.M.T.  >, 
evidence  of  a  limited  isallobaric  contribution  was  given  by  the 
geostrophic  crossfrontal  components  in  better  agreement  with  the 
observed  frontal  velocity  than  the  actual  wind  crossfrontal 
components. 

3.2  15  November  1988  Case  .Study 

At  1200  G. M.T.  an  intense  shortwave  moved  out  of  the  southwestern 
United  States  (note  the  strength  of  the  500  mb  height  falls  over 
northern  New  Mexico,  Fig.  39)  into  a  barocllnlc  environment  favorable 
for  rapid  surface  cyclogenesis  (Fig.  40).  This  well  defined 
barocllnlc  zone  over  Kansas  and  Nebraska  was  predominately  due  to 
vigorous  warm  air  adveetlon  In  the  lower  levels  (Fig.  40).  This  warm 
air  adveetlon  over-  the  past  12  hours  moved  the  warm  front  from  Kansas 
to  the  1600  G.M.T.  position  over  Minnesota  and  Wisconsin  (Fig.  41). 
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As  indicated  by  the  NMC  1200  G.M.T.  850  mb  chart  (Fig.  40),  cold  air 
advection  extended  from  South  Dakota  to  western  Nebraska.  By  the 
beginning  of  the  case  study  <1600  G.M.T. ,  Fig.  41),  cold  air 
advection  promoted  the  development  of  the  cold  front  from  the 
Oklahoma  Panhandle  to  northern  New  Mexico  and  was  partially 
responsible  for  the  cold  front  from  Kansas  to  eastern  Nebraska.  The 
upper  level  shortwave  concurrent  with  the  baroclinic  zone  facilitated 
the  development  of  an  intense  surface  cyclone  over  southwestern 
Kansas  that  moved  into  northeast  Kansas  by  2300  G.M.T.  (Fig-  41). 

The  cold  front  extending  from  Oklahoma  into  Texas  at  1600  G.M.T. 

(Fig.  41)  stemmed  from  advection  of  dry  air  <Fig.  40)  and  subsequent 
frontal  development  due  to  an  abrupt  change  in  moisture. 

In  addition  to  surface  cyclone  development,  this  case  study 
differs  from  the  31  January  case  study  in  two  important  respects: 
First,  widespread  convective  precipitation  was  indicated  throughout 
the  entire  period  of  consideration  <Figs,  42,  43,  and  44).  Many 

tornadoes  were  reported  in  the  impressive  line  of  thunderstorms  that 
materialized  from  Iowa  to  Texas.  Thundersnows  were  reported  in  the 
isolated  convection  over  central  and  northwestern  Kansas.  Second, 
rapidly  changing  D  value  tendencies  at  the  surface  created  a  myriad 
of  isallobaric  forcing  patterns,  while  the  :3L  January  case  study 
exhibited  a  more  consistent  isallobaric  pattern.  In  this  case  study, 
cold  air  advection  initially  resulted  in  D  value  rises  (1600  to  1800 
G. M.T. )  behind  the  cold  front  extending  from  Kansas  to  Nebraska.  As 
the  surface  cyclone  moved  northeast,  out  of  southwestern  Kansas,  the 
isallobaric  acceleration  on  the  actual  winds  in  the  vicinity  of  this 
front  was  in  response  to  l)  value  falls.  For  example,  a  pressure  fall 
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Ftg.  42.  National  Meteorological  Center  Radar  Summary  at  1735  G.M.T. 

Hatched  areas  are  precipitation.  Contours  are  at  intensity 
levels  one,  three,  and  five.  Area  movements  (kts)  are 
given  by  arrows  with  pennants.  Cell  movements  <kts>  are 
given  by  arrows  with  numbers.  Precipitation  bases  and  tops 
( X  10*  feet)  are  overlined  and  underlined,  respectively. 
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Fig.  44.  Same  as  Fig.  42,  except  at  2335  G.M.T.  . 


57 

of  3.6  mb  was  reported  at.  Topeka,  Kansas  from  1700  to  1800  G.  M.T. 

Given  the  movement  and  continued  development  of  the  surface  cyclone 
over  Kansas,  the  cold  front  extending  from  Kansas  to  Nebraska 
displayed  an  anomalous  movement.  I sal  lobar ic  forcing  appears  to 
provide  a  partial  explanation  for  the  movement  of  this  front. 

In  response  to  the  D  value  tendencies  over  Nebraska  and  Kansas 
(Fig.  45),  the  1600  G.M.T.  agoostrophic  winds  (Fig.  46)  and  frontal 
velocities  (Fig.  47)  hint  of  an  isallobaric  contribution  over  frontal 
section  G- 1 .  Despite  the  presence  of  strong  isallobaric  crosofrontal 
components,  an  isallobaric  influence  on  the  actual  wind  crossfrontal 
components  appears  limited  to  frontal  segment  H.  Fr ontogone t ica 1 
forcing  due  to  the  isallobaric  wind  <Fig.  48)  indicated  a  maxima  of 
f ran toge nes is  concurrent  with  the  I)  value  rise  center  (Fig.  45> 
behind  frontal  section  G-K.  A  similar  f rontogenes is  pattern 
displayed  by  the  f rontogenet i ca 1  forcing  due  to  the  actual  wind  (Fig. 
49)  alludes  to  isallobaric  forcing  as  a  contributor. 

The  actual  wind  f rontogenet ica 1  forcing  would  contribute  to 

slower  movement  of  frontal  section  G-K,  but  the  magnitude  of 

o  **/  %/ 

temperature  gradient  realignment  is  small  (1  K  hr  100  km  ).  While 
the  actual  wind  crossfrontal  components  provide  a  more  consistent  (in 
the  direction  of  motion)  indication  of  movement  for  frontal  section 
G  I,  an  isallobaric  con  t.  r  1  but  i  on  is  difficult  to  diagnose  at  frontal 
segments  G  and  I.  The  limited  area  of  isallobaric  influence  on  the 
actual  crossfrontal  components  is  demonstrated  by  geostrophic 
crossfrontal  components  in  better  agreement  with  the  observed  frontal 
motion  for  frontal  section  A-F.  in  contrast,  the  actual  wind 
crossf r onta 1  components 


conform  better  to  t  lie  observed  frontal 
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Frontal  velocity  map  for  the  period  1600-1700  0. M. T.  15 
November  1988.  Frontal  velocities  <m  s’*)  are  indicated  by 
the  larger  numbers  and  correspond  to  the  frontal  segments 
Indicated  in  parenthesis  by  the  frontal  velocity.  Arrows 
indicate  the  approximate  direction  of  motion.  Small  numbers 
indicate  the  crossfrontal  component  (ra  s’*)  by  the 
actua 1/geostroph ic/isa 1  lobar ic  winds.  Positive  (negative) 
cross  frontal  components  indicate  flow  in  the  same 
(opposite)  direction  of  frontal  motion. 


60 


Fig.  48.  Frantagenet ica 1  farcing  due  ta  the  Isa  1 lobar lc  wind  at 
1600  G.M.T,  15  November  1988  (X  103®K  kmf hr) .  Positive 
(negative)  values  represent  f rontogenesis  (f rontolysis)  . 
Surface  frontal  positions  are  superimposed. 
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velocity  than  the  geostrophic  crossfrontal  components  over  frontal 
section  J-T  were  isallobaric  contributions  are  nonexistent. 

Continued  development  of  the  surface  cyclone,  as  it  moved  into 
north-central  Kansas,  facilitated  a  different  D  value  tendency  field 
(.Fig.  50)  by  1800  G.M.T.  The  1800  G.M.T.  ageostrophic  winds  (Fig. 

51)  over  frontal  section  H-L  in  conjunction  with  the  magnitude  of  the 
acttial,  geostrophic,  and  isallobaric  crossfrontal  components  <Fig. 

52)  initially  allude  to  an  isallobaric  contribution,  but  a  frictional 
contribution  is  suggested  as  well.  Advocation  of  isallobaric 
influence  on  the  actual  wind  crossfrontal  components  at  frontal 
segments  J,  K,  and  L  is  supported  through  similar  magnitudes,  but  in 
the  opposite  direction  of  actual  versus  geostrophic  crossfrontal 
components.  The  same  comparison  on  frontal  segments  H  and  I  is 
ambiguous  due  to  smaller  (in  apposite  direction)  actual  versus 
geostrophic  crossfrontal  components.  Similar  f rontogenet ical  forcing 
by  the  1800  G.M.T.  actual  and  isallobaric  winds  (Figs.  53  and  54, 
respectively)  illustrated  patterns  favorable  for  "apparent" 
propagation  over  frontal  section  II—  J  .  As  opposed  to  the  31  January 
case  study,  the  isallobaric  contribution  to  "apparent”  propagation 
included  both  f rontogenesis  ahead  of  the  front  as  well  as  the 
frontolysis  behind  the  front.  S ign i f ica nt 1 y ,  the  isallobaric  forcing 
occurred  in  a  D  tendency  pattern  where  falls  ahead  and  behind  the 
front  are  observed.  Recall  that  in  the  31  January  case  study,  the 
suggested  isallobaric  contribution  to  "apparent"  propagation  was 
limited  to  the  frontolysis  behind  the  front. 

As  observed  for  1600  G. M.T. ,  the  geostrophic  crossfrontal 
components  are  not  supportive  of  the  observed  motion  for  section 


Fig.  54.  Same  as  Fig.  48,  except  at  1800  G.M.T. . 
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H-K.  Unlike  1600  G.M.  T.  where  lsallobaric  influence  was  limited  to 
one  actual  wind  crossf rontal  component,  both  the  1800  G.M. T.  actual 
crossfrontal  components  and  f rontogenet ica 1  forcing  are  supportive  of 
the  observed  frontal  motion.  Actual  wind  crossfrontal  components  in 
better  agreement  with  the  observed  frontal  motion  are  not  limited  to 
frontal  section  H-K,  but  it  is  over  this  frontal  section  that  an 
lsallobaric  influence  is  suggested. 

The  appearance  of  a  D  value  rise  rise  maxima  in  the  2000  G.M.T.  D 
value  tendencies  introduced  a  new  Lsallobaric  forcing  field  (Fig.  55) 
behind  frontal  section  l-L  (Fig.  56).  lsallobaric  contributions  to 
actual  crossfrontal  components  (Fig.  57)  a  e  suggested  over  frontal 
section  J-L.  Frontogene t J ca i  forcing  due  to  the  lsallobaric  wind 
(Fig.  58)  now  demonstrated  f rontogones is  behind  and  f rontolysis  ahead 
of  frontal  section  J-K,  while  frontal  segments  L  and  M  were  still 
indicative  of  "apparent"  propaga t ion .  The  f rontogenet ica 1  forcing 
due  to  the  actual  wind  <Flg.  56)  also  demonstrated  a  similar  pattern, 
but  suggests  that  the  amount  of  f rontogenet ica 1 ly  induced  motion 
decreased  across  frontal  section  J-K.  In  summary,  the  appearance  of 
I)  value  rises  behind  frontal  section  I-K  decreased  the 

f rontogenet ica 1  forcing  favorable  for  "apparent"  propagation.  Where 
the  pattern  of  D  value  falls  remained  across  life  front,  "apparent" 
propagation  was  still  evident.  While  decreased  amount  of 
f ron togenet ica 1 1 y  induced  motion  across  frontal  section  I-K  raises 
doubts  that  the  observed  frontal  motion  is  well  diagnosed  by  the 
actual  wind  crossf ron ha  I  components  and  1 r on togonn t i ca 1  forcing,  the 
combination  provides  a  bother  indication  ol  the  observed  frontal 
motion  than  the  geostrnphie  crossfrontal  components  alone. 


In  this  case  study,  the  anomalous  movement  displayed  by  the  front 
over  eastern  Nebraska  was  fairly  well  diagnosed  by  the  actual  wind 
crossfrontal  components  and  f rontogenet  ica 1  farcing.  The  isallobaric 
f rontogenet ica 1  forcing  as  a  result  of  D  value  rises  behind  the  front 
and  D  value  falls  ahead  of  the  front  for  both  1600  and  2000  G. M. T. 
was  not  favorable  for  "apparent"  propagation.  Similarly  observed 
patterns  of  f ron togenet ica 1  forcing  duo  to  the  actual  wind  allude  to 
isallobaric  contributions,  but  the  results  of  the  15  November  and  31 
January  case  studies  are  not  supportive  of  "apparent”  propagation  due 
solely  to  isallobaric  forcing  in  the  hypothesized  D  tendency  field 
shown  in  Fig.  1.  In  contrast,  isa 1  lobar ica 1 1 y  influenced 
f rontogenet ica 1  f^rcing  favorable  for  "apparent"  propagation  appeared 
when  larger  D  value  falls  were  indicated  ahead  of  the  front. 
Unfortunately,  the  Isallobaric  influenced  f rontogenet i ca 1  forcing 
occurred  during  the  period  when  thunderstorms  were  present  in  the 
area.  The  effects  of  the  convective  precipitation  over  the  region 
introduce  another  variable  to  the  f rontogene t ica 1  forcing  equation  by 
virtue  of  changed  temperature  gradients  not  due  to  the  winds  alone. 
Due  to  a  somewhat  consistent  isallobaric  influence  on  the  actual  wind 
crossfrontal  components  for  the  front  over  Nebraska  and  Iowa,  the 
isallobaric  contribution  to  surface  frontal  motion  appears  more 
significant  than  in  the  31  January  case  study.  The  question  that 
immediately  comes  to  mind  is  why  did  Isallobaric  influences  play  a 
more  significant  role  for  this  front  as  opposed  to  the  other  fronts 
in  this  case  study?  In  addition,  why  didn't  the  31  January  case 
study  with  substantially  greater  I>  value  tendencies  show  a  stronger 
level  of  Influence?  The  answer  may  come  from  the  other  possible 


Fig.  55.  Same  as  Fig 


except  at  2000  G. M. T 
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contributors  to  surface  frontal  motion,  such  as  the  crossfrontal 
components  of  the  upper  level  winds.  In  the  15  November  case  study, 
the  upper-  level  winds  (Fig.  60)  appear  parallel  to  the  frontal 
surface.  If  the  upper  level  winds  were  parallel  to  the  frontal 
surface,  the  surface  frontal  mot;  ion  would  be  controlled  to  a  larger 
extent  by  surface  forces  than  would  normally  be  the  case.  This  is 
only  a  theory  and  without  more  in-depth  analysis  cannot  be 
subs t ant  fated . 

3.3  7  January  1989  Case  Study 

The  7  January  case  study  provides  another  analysis  of  isal lobaric 
effects  on  surface  frontal  motion  associated  with  a  developing 
surface  cyclone.  This  case  study  exhibits  two  primary  synoptic 
differences  from  the  previous  case  study:  First,  the  combination  of 

slow  movement  by  the  developing  extratropical  cyclone  and 
availability  of  cold  air  (Fig.  61).  Second,  a  longer  wavelength  by 
the  upper  level  shortwave  (Fig.  62)  conducive  to  benign  las  compared 
to  the  15  November  case  study)  surface  cyclone  development  through 
baroclinic  instability.  Subdued  surface  cyclogenesis  concurrent  with 
an  ample  supply  of  cold  air  produced  generally  constant  isal lobaric 
forcing  (similar  to  the  31  January  case  study)  on  the  surface 
fronts.  D  value  falls  and  rises  were  observed  ahead  and  behind  of 
the  eastern  cold  front,  respectively  extending  from  Minnesota  to 
Texas  at  1.700  G.M.T.  (Fig,  63).  in  contrast,  both  cold  fronts  had 
similar  origins  as  the  cold  fronts  in  the  15  November  case  study. 
Initially,  the  eastern  cold  front  was  the  result  of  an  abrupt  change 
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in  moisture.  The  western  cold  front  at  L700  G.M.T.  (Minnesota  to  New 
Mexico,  Fig.  64),  originated  from  the  cold  air  advection  behind  the 
surface  cyclone. 

Precipitation,  limited  at  the  inception  of  the  case  study  (Fig. 
65),  rapidly  developed  (Figs.  66  and  67)  and  by  2235  G.M.T.  (Fig.  68) 
substantial  areal  coverage  by  showers  and  thunderstorms  was 
indicated.  Stratiform  precipitation  was  widespread,  but  small  areal 
coverage  on  the  NMC  Radar  Summaries  can  be  attributed  to  the  fact 
that  it  was  snow.  Initially  therefore,  it  is  passible  to  analyze  the 
isallobaric  forcing  on  the  surface  frontal  motion  and  remove  the 
unknown  effects  of  convective  precipitation  on  f rontogenet ica 1 
fore i ng . 

The  1800  G.M.T.  D  value  tendencies  (Fig.  69)  demonstrate  the 
hypothesized  pattern  favorable  for  an  isallobaric  contribution  to 
frontal  motion  over  cold  frontal  sections  A- 1  and  a-d  (Fig.  70). 

Even  though  the  largest  isallobaric  crossfrontal  components  (Fig.  71) 
are  Indicated  on  frontal  section  E~J ,  Isallobaric  influence  on  the 
actual  crossfrontal.  components  is  limited  4o  frontal  segments  a,  b, 
and  d.  Frontoge net ica 1  forcing  by  the  1800  G.M.T.  actual  and 

isallobaric  winds  (Figs.  72  and  73,  respectively)  suggest  "apparent” 
propagation  by  isallobaric  forcing  over  frontal  sections  a-c,  f-h, 
and  to  a  lesser  extent  slowed  movement  over  frontal  section  E-J. 
Combination  of  the  f ron t oge ne t ica 1  forcing  due  to  the  actual  wind  and 
the  actual  wind  crossfrontal  components  produces  a  limited  area  where 
a  better  diagnosis  of  surface  frontal  motion  than  the  geostrophic 
crossfrontal  components  (segments  a,  b,  and  d)  is  achieved.  However, 
these  are  the  segments  where  an  isallobaric  influence  is  possible. 
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Fig.  70.  1800  G.M.T.  ageostrophlc  winds  plotted  at  the  grid  points 

used  to  determine  the  crossfrontal  components  by  the 
actual,  geostrophic,  and  lsallobarlc  winds.  Letters 
designate  frontal  segments  on  superimposed  fronts. 
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Fig.  71.  Frontal  velocity  map  for  the  period  1800-1900  G.M.T.  7 

January  1989.  Frontal  velocities  (m  s‘)  are  Indicated  by 
the  larger  numbers  and  correspond  to  the  frontal  segments 
indicated  In  parenthesis  by  the  frontal  velocity.  Arrows 
indicate  the  approximate  direction  of  motion.  Small  numbers 
Indicate  the  crossfrontal  component  <m  s”*>  by  the 
aetua 1/geostrophic/isa 1  lobar ic  winds.  Positive  (negative) 
cross  frontal  components  indicate  flow  in  the  same 
(opposite)  direction  of  frontal  motion. 


Frontogenet i ca 1  forcing  due  to  the  actual  wind  at 
1800  G.M.T.  7  January  1989  v X  103°K  Km 1  hr'*  >  .  Positive 
(negative)  values  represent  f rontogenesis  (frontolysis) 
Surface  frontal  positions  are  superimposed. 


Fig.  74.  Sa me  as  Fig 


except  at  1900  G.M.T 
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The  actual  wind  f rontogenet ica 1  forcing  is  cont r i but i ng  to  an 
eastward  propagation  of  frontal  section  i-h,  but  the  suggested 
isal lobar ic  influence  on  the  actual  wind  f rontogenet ica 1  forcing 
promoting  the  propagation  occurs  in  a  gradient  of  D  value  rises.  In 
contrast,  the  "apparent"  propagation  potentially  due  to  isal lobaric 
influence  on  the  actual  wind  f rontogenet leal  forcing  over  frontal 
section  a-c  lies  in  I)  tendencies  hypothesized  as  favorable  for  this 
phenomenon.  Similar  results  are  observed  for  the  1900  G.M. T. 
lsallobaric  (Fig.  74)  and  actual  wind  (Fig.  75)  f  rontogenet  ical 
forcing  with  only  minor  changes  observed  in  the  D  tendency  field 
(Fig.  76).  These  differences  in  f rontogenet ica 1  forcing  versus  the  D 
tendencies  were  due  to  the  temperature  tendencies  term.  Notice  the 
difference  in  the  1900  G.M. T.  gradient  of  D  tendencies  across  the 
frontal  section  over  Minnesota  and  the  frontal  section  extending  from 
central  Iowa  to  northern  Missouri  (Fig.  76).  The  1900  G.M.T. 
temperatures  fitted  to  the  grid  (Fig.  77),  illustrate  that  the 
stronger  gradient  of  D  change  values  (and  therefore  isal lobaric 
winds,  Fig.  76)  occurs  in  the  maximum  temperature  gradient.  As 
demonstrated  in  the  31  January  case  study,  this  combination  promotes 
f rontoge nes i s  due  to  gradient  of  advection  portion  of  the  temperature 
tendencies  terra.  In  contrast,  over  southwestern  Missouri  the  maximum 
temperature  gradient  (Fig.  77)  is  over  the  D  rise  maxima  and  a  broad 
area  ol  tici.tolysis  occurs.  The  weaker  temperature  gradient  combined 
with  a  weaker  D  tendency  gradient  associated  with  the  front  over 
Minnesota  allowed  a  dominant  divergence  terra  and  subsequent 
"apparent"  propagation  was  indicated  as  in  the  hypothesized  D 
tendency  field  (Fig.  1). 


Fig.  75.  Same  as  Fig.  72,  except  at  1900  O.M.T 


Fig.  76.  Same  as  Fig.  69,  except  at  1900  G.M.T 
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Fig.  77.  Temparatures  ,0K)  fitted  to  the  grid  points  for  1900  G.M.T 
7  January  1989.  Dashed  lines  are  isotherms  analysed  at  4® K 
intervals.  Surface  frontal  postions  are  superimposed. 


Fig.  78. 


Same  as  Fig.  76,  except  at  2000  G.  M.  T.  . 
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By  20  00  G.M.  T.  ,  the  D  value  rise  gradient  (Fig.  78)  now  extended 
across  frontal  section  B-K  (Fig.  79).  Isallobaric  influence  on  the 
actual  wind  crossfrontal  components  is  suggested  on  frontal  segments 
a-d,  1-n,  and  H,  J,  K,  and  L  (Figs.  79  and  80).  Notice  that  this 
hour  had  the  rare  occurrence  of  an  isallobaric  influence  where  the 
maximum  isallobaric  crossfrontal  component  was  indicated  (segment  H, 
Fig.  79).  In  addition,  the  first  occurrence  of  an  isallobaric 
influence  on  the  actual  wind  crossfrontal  component  associated  with 
the  warm  front  is  suggested.  Frontogenet ica 1  forcing  due  to  the 
isallobaric  wind  (Fig.  81)  demonstrated  a  change  similarly  observed 
in  the  31  January  case  study.  In  particular,  the  sudden  appearance 
of  a  large  area  of  frontolysis.  While  frontolysis  was  occurring 
across  both  of  the  cold  fronts  in  previous  hours,  its  location  was 
farther  south  and  did  not  encompass  the  area  indicated  at  2000 
G.M.T.  Recall  In  the  31  January  case  study,  that  the  frontolysis 
was  due  to  a  repositioned  D  value  rise  center  over  the  maximum 
temperature  gradient  and  subsequent  dominance  by  the  divergence 
term.  In  this  case,  the  reason  was  due  to  a  combination  of  terms. 

In  addition,  the  decrease  In  the  2000  G.M.T.  I)  tendency  gradient 
(Fig.  78)  over  the  location  of  the  maximum  temperature  gradient, 
provides  a  partial  explanation  for  the  absence  ol  f rontogenesis.  It 
is  interesting  to  note  the  decreased  magnitude  of  D  value  tendencies 
and  a  decreased  isallobaric  contribution  <as  compared  to  the  31 
January  case  study)  to  the  actual  wind  f ron togenet ica 1  forcing  (Fig. 
82).  Isallobaric  contributions  to  frontal  motion  are  again 
suggested,  but  the  geostrophic  wind  crossfrontal  components  are  more 
useful  for  diagnosing  most  of  the  observed  frontal  motion  from  2000 


Fig.  70.  Same  as  Fig.  70,  except  at  2000  G.M.T. 


Fig.  80.  Same  as  Fig,  71,  except  at  2000  G.M.T. 
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G.  M.  T.  to  2100  G.  M.  T. 

This  case  study  was  similar  to  the  two  previous  case  studies  in 
that  limited  and  somewhat  random  areas  of  isallobaric  influences  on 
the  actual  wind  crossfrontal  components  were  indicated.  Similarities 
to  the  15  November  case  study  include  the  following:  First,  only  one 

section  of  a  cold  front  displayed  a  consistent  isallobaric  influence 
on  the  actual  wind  crossfrontal  components.  Furthermore,  the 
observed  movement  of  this  frontal  section  was  in  better  agreement 
with  actual  than  geostrophic  crossfrontal  components.  Second, 
although  the  D  tendency  field  was  different,  the  frontal  section 
where  an  isallobaric  influence  on  the  actual  wind  crossfrontal 
components  was  suggested  also  indicated  an  isallobaric 
frontogenetical  farcing  contribution  which  promoted  ’’apparent” 
propagation.  The  most  significant  similarity  to  the  31  January  case 
study  was  the  appearance  of  isallobaric  induced  frontolysis  in  the 
frontogenetical  forcing  due  to  the  actual  wind.  In  both  case 
studies,  this  was  due  to  the  D  value  rise  center  positioned  over  the 
maximum  temperature  gradient.  However,  the  decreased  magnitude  of 
Isallobaric  accelerations  in  this  case  study  still  contributed  to 
similar  patterns  in  the  actual  wind  frontogenetical  forcing.  It  is 
passible  therefore,  that  the  amount  of  Isallobaric  influence  on  the 
the  actual  wind  frontogenetical  forcing  is  directly  related  to  the 
strength  of  the  isallobaric  acceleration.  This  hypothesis  has  two 
assumptions:  Firstly,  that  the  magnitude  of  the  temperature 

gradients  are  similar.  Secondly,  other  terms  comprising  the 
ageostrcphic  wind  are  similar  in  magnitude.  While  the  first 
assumption  is  fulfilled  in  this  case,  the  second  assumption  is  by  no 
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means  a  certainty.  However,  it  is  possible  that  friction  at  the 
surface  is  so  dominant  that  the  magnitudes  of  the  other  accelerations 
(besides  the  isallobaric  acceleration)  comprising  the  ageostrophic 
wind  are  insignificant. 

4.  Summary  and  Suggestions  For  Future  Research 

Despite  the  fact  that  the  isallobaric  influence  of  the  actual 
wind  crossfrontal  components  had  f rontogenet ica 1  forcing  had  to  be 
subjectively  diagnosed,  this  research  has  provided  observational 
evidence  that  isallobaric  forcing  can  be  a  contributor  to  surface 
frontal  motion  and  f rontogenesis .  Indications  of  an  isallobaric 
contribution  to  surface  frontal  motion  through  comparison  of  actual 
and  geostrophic  crossfrontal  components  are  constrained  by  three 
assumpt ions : 

1)  The  assumption  of  a  front  as  a  material  surface. 

This  assumption  has  been  debated  since  the  inception  of  fronts 
and  to  date  neither  frontal  theory  nor  observation  has  provided  an 
answer.  As  noted  by  Reeder  and  Smith  ' 1988),  the  assumption  that  a 
front  is  a  material  surface  is  useless  in  the  case  of  a  frontal  zone 
which  moves  faster  than  the  measured  wind  ahead  or  behind  the  front 
at  any  le ve  1 . 

2)  The  assumption  that  the  winds  at  the  grid  points  used  to 
compute  the  crossfrontal  components  are  representative  of  the  force 
exerted  by  the  wind  on  the  front. 

In  particular,  is  the  force  imparted  by  the  wind  accurately 
represented  on  the  frontal  scale.  Unfortunately,  the  answer  is 
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la  probably  no.  The  Inability  to  determine  the  significance  of 
crossfrontal  components  on  frontal  motion  opens  the  question  of 
frontal  location  versus  the  location  of  the  wind  that  gives  a  true 
measure  of  the  force  exerted  on  the  front.  However,  Brundidge's 
(1965)  results  suggest  that  the  most  accurate  measure  of  the  force 
imparted  on  a  front  by  the  wind  is  close  behind  the  frontal  surface 
and  the  grid  in  this  study  does  not  have  the  resolution  needed  to 
represent  the  winds  on  the  frontal  scale. 

3)  The  assumption  that  the  frictionless  isallobaric  wind  is 
representative  of  the  isallobaric  wind  at  the  surface. 

To  some  degree  this  assumption  is  valid;  however,  in  theory  the 
amount  of  difference  between  the  two  can  be  substantial  (Young, 

1973).  The  fact  that  stronger  actual  than  geostrophic  wind 
crossfrontal  components  were  probably  due  to  the  isallobaric 
acceleration  does  suggest  that  there  can  be  an  isallobaric 
contribution  to  frontal  motion.  However,  the  limited  areas  of 
isallobaric  Influence  on  the  actual  wind  crossfrontal  components 
suggests  that  the  isallobaric  contribution  to  frontal  motion  is 
1 i mi ted . 

The  actual  wind  f rontogenet leal  forcing  results  of  this  study 
have  demonstrated  the  occurrence  of  f r ontogene t  ica  1  forcing  patterns 
favorable  for  "apparent"  propagation.  Furthermore,  this  study  has 
subjectively  shown  that  isallobaric  forcing  can  be  a  significant 
contributor  to  "apparent"  propagation,  but  the  diagnosed  isallobaric 
influence  on  the  actual  wind  f rontogenet ica 1  forcing  is  constrained 
by  the  third  assumption  above  as  well.  In  the  15  November  case 
study,  convective  precipitation  could  have  contributed  to  the 
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"apparent"  propagation  demonstrated  by  the  actual  wind 
f rontogenet ica 1  forcing.  This  study  has  shown  that  the  temperature 
gradients  associated  with  a  front  can  go  through  a  process  of 
redevelopment.  Obviously,  three  case  studies  do  not  give  a  full 
measure  of  the  isallobaric  influence  on  the  actual  wind 
f rontogenet ical  forcing,  but  the  results  do  allude  to  pressure 
tendency  patterns  most  favorable  for  "apparent”  propagation  in  the 
actual  wind.  Provided  that  the  level  of  Isallobaric  influence  is  the 
same  as  observed  in  this  study,  these  are  as  follows:  First,  large 

pressure  falls  ahead  of  the  maximum  temperature  gradient  with  smaller 
pressure  falls  on  or  behind.  Second,  a  weak  gradient  of  pressure 
rises  behind  and  pressure  falls  ahead  of  the  maximum  temperature 
gradient,  respectively.  Finally,  it  is  possible,  that  when  a 
significant  pressure  rise  center  moves  over  the  maximum  temperature 
gradient  frontolysls  will  be  indicated  by  the  actual  wind 
f  rontogene  t  ica  1  forcing.  Whether  or  not  the  froiitogenes  is  ahead  of 
the  temperature  gradient  will  occur  is  uncertain. 

Given  the  high  degree  of  uncertainty  associated  with 
crossfrontal  components,  the  f rontogenet ica 1  forcing  results  provide 
the  most  promising  area  for  future  research.  Recall  that  the 
f rontogenet ica 1  forcing  computations  in  this  study  considered  only 
temperature.  Since  moisture  amounts  can  affect  the  density  of  two 
different  air  masses,  a  follow-on  study  which  includes  the  effects  of 
moisture  by  use  of  the  virtual  temperature  could  easily  be  done.  A 
less  subjective  method  to  define  the  Isallobaric  influence  on  the 
actual  surface  wind  f rontogenet lea 1  forcing  would  be  to  Include  the 
effects  of  friction  on  the  Isallobaric  wind.  In  addition,  a  wider 
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range  of  pressure  tendency  fields  should  be  covered  to  determine  the 
limitations  of  the  isallobaric  influence  on  the  actual  wind  versus 
the  magnitude  of  the  pressure  tendency.  Young  (1973)  has  established 
the  theoretical  groundwork  for  inclusion  of  the  frictional 
isallobaric  wind.  However,  a  study  of  this  nature  would  still  have 
uncertainty  due  to  the  other  terms  that  comprise  the  ageostrophic 
wind.  Finally,  a  study  to  determine  the  velocity  induced  on  a  front 
by  the  f rontogenet ica 1  forcing  process  would  be  constrained  by  the 
uncertainty  of  temperature  gradient  velocity  versus  frontal  velocity, 
but  could  establish  a  range  of  velocities  versus  the  magnitudes  of 
the  temperature  gradient. 

APPENDTX  A  -  Data  Checking  Methods  Used 

Screening  for  errors  in  the  data  set  consisted  of  a  three  step 
process.  First,  all  station  data  was  compared  to  predetermined 
boundary  values  for  each  element.  The  high  and  low  boundary  values 
used  we~e  as  follows:  75  and  U  *  F  tor  temperature,  60  and  0°  F  for 
dewpoint,  360  and  0  degrees  for  wind  direction,  25  and  0  kts  for  wind 
speed,  and  30.50  and  29.80  inches  for  altimeter  settings.  Values 
outside  the  boundary  conditions  were  identified  and  checked  against 
the  original  and  corrected  airway  reports.  Second,  temperatures, 
dewpoints,  and  altimeter  settings  at  every  observation  location  were 
compared  to  the  same  elements  at  all  other  reporting  stations  within 
a  130  km  radius  from  the  observation  location.  Values  at  the 
observation  location  outside  the  specified  range  for  each  element 
were  identified  and  checked  against  the  original  teletype  reports. 
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The  ranges  used  were  as  follows:  5®  F  for  temperature,  5  °F  for 
dewpoint,  and  .05  inches  for  altimeter  set. ting.  Due  to  the  high 
variability  in  wind  direction  and  speed  by  local  effects  alone,  wind 
direction  and  speed  were  not  checked  in  this  manner.  Finally, 
temperature,  dewpoint,  wind  direction,  wind  speed,  and  altimeter 
setting  values  at  each  observation  location  were  identified  and 
checked  if  the  difference  was  greater  than  '5  0  F ,  4®  F,  35°,  5  kts, 

and  .05  inches,  respectively  from  the  preceding  hour. 

A  more  subjective  correction  of  improperly  dissmenated  elements 
was  required  in  same  cases.  Obvious  dissemination  errors  were  easily 
identified  and  corrected.  Far  example,  the  altimeter  setting 
reported  as  28.67"  in  a  Intermediate  hour  and  reported  as  29.64”  in 
the  preceding  hour  and  29.70"  in  the  third  hour.  More  difficult  to 
identify  as  errors  were  the  not  impossibly  large  changes  of  a 
particular  element.  For  example,  large  pressure  changes  were  checked 
against  required  remarks,  pressure  values  for  the  hours  on  either 
side,  and  against  pressure  changes  by  nearby  locations.  Overwhelming 
evidence  was  required  to  remove  data  that  fell  into  this  category. 
These  data  checking  methods  do  not  guarantee  an  error  free  data  set 
as  smaller  errors  by  dissemination  and  data  entry  could  pass  through 
the  error  checking  routines  and  introduce  noise  into  the  results. 

In  some  case  studies,  the  strengtn  of  the  computed  isaliobaric 
wind  raised  doubts  about  the  validity  of  the  method.  Confirmation  of 
the  isaliobaric  wind  values  computed  was  achieved  by  another  computer 
program  using  the  change  of  pressure  formula  <, Fq .  17)  for  the 

isaliobaric  wind.  Perfect  agreement  was  found  at  most  grid  points. 
Small  differences  in  the  wind  direction  and  wind  speed  did  occur  at 
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same  grid  points,  These  differences  were  mast  likely  due  to  rounding 
by  the  computer  as  the  the  wind  direction  at  each  grid  point  was 
rounded  to  the  nearest  10®.  The  largest  differences  noted  were  on 
the  order  of  30®  and  4  kts  far  the  wind  direction  and  wind  speed, 
respectively.  Directional  differences  of  this  magnitude  occurred 
where  the  isal lobar ic  wind  magnitudes  were  weak  and  speed  differences 
occured  where  the  magnitudes  were  exceptionally  large. 

APPENDIX  B  -  Description  of  the  Cressman  <1959)  Weighting  Scheme 

Since  an  objective  analysis  produces  an  approximation  of  the 
value  of  an  observed  or  computed  parameter  at  grid  points, 
consideration  of  the  error  produced  by  the  objective  analysis  scheme 
is  necessary.  The  computer  code  used  for  the  Cressman  (1959) 
weighting  scheme  is  a  modified  form  of  the  code  presented  by  Inman 
<1970).  A  general  discussion  will  be  used  to  qualitatively  asses  the 
error  introduced  by  the  objective  analysis  method.  Let  T*  represent 

the  temperature  of  a  given  observation  location  where  s  is  the  scan 

-  f .  . 

number.  The  correction  applied  to  the  grid  point  i,j  under 

consideration  by  the  temperature  of  the  observation  location  is 
calculated  by 


cB  .  -  W  Tb/V. 

<• .  J 


<B1> 


V  in  Eq.  <B1)  is  a  distance  dependent  weight  factor  that  is 
determined  by  the  distance  of  the  observation  location  from  the  grid 
point  using  the  equation  below. 
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V  =  (CR2  -  d2)CR2  +  d2>  4>, 


(B2) 


where  d  is  the  distance  between  the  grid  point  and  the  observation 
location  and  R  is  the  radius  of  influence  specified.  The  magnitude 
of  the  weighting  factor  V  from  Eq.  (B2)  based  on  the  distance  between 
the  grid  point  and  the  observation  location  is  depicted  in  Fig.  Bl. 
The  total  correction  applied  to  the  grid  point  is  a  weighted  average 
of  all  the  observation  locations  within  the  radius  of  influence. 

When  the  distance  <d)  between  the  grid  point  and  the  observation 
location  is  larger  than  R,  the  station  element  will  not  influence  the 
correction  applied  to  the  grid  point  <W  =  0).  When  the  correction 
due  to  all  stations  within  the  radius  of  influence  has  been  computed, 
the  estimate  of  temperature  at  the  grid  point  is  calculated  by 


T® 

>- .  j 


__  ^.s  -  1 


L .  J 


+  C 


>-  »  J 


<B3 ) 


For  the  first  scan,  T  represents  the  first  guess  field.  As 
McDonnell  <1967)  noted,  first  guess  fields  are  not  absolutely 
necessary  for  a  successful  application  of  the  technique  where 
sufficient  data  are  available.  Therefore,  a  zero  first  guess  field 
was  used.  Before  the  next  scan  is  executed,  the  current  estimate  of 
T is  made  at  the  observation  location.  The  estimate  is  made  from 
values  at  the  four  surrounding  gridpoints  using  the  bilinear 
interpolation  formula  as  given  by  Inman  <1970), 


T8  =  T®  +  (T®  -  T® >  Ay/b  +  (T®  -  T®)Ax/a  - 

pl  4  1  J  2  1 

CT2  -  T®  +  T®  -  T® ) AxAy/ab, 


<B4> 
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where  Is  the  estimate  of  temperature  at  the  observation  location, 
T*  with  a  numbered  subscript  is  the  temperature  at  one  of  the  four 
surrounding  grid  paints,  Ax  is  the  distance  in  the  x-direction  from 
the  y-coordinate  of  grid  points  one  and  four,  Ay  is  the  distance  in 
the  y-direction  from  the  x-coordinate  of  the  grid  paints,  a  is  the 
distance  between  grid  points  along  the  x-axis,  and  b  is  the  distance 
between  grid  paints  along  the  y-axis.  Illustration  of  the  terms  in 
Eq.  <B4>,  presented  in  Fig.  B2 ,  will  help  clarify  the  subsequent 
discussion.  The  formula  assumes  a  square  grid  is  used.  For  the 
lat itude/longitude  grid  used,  the  y-direction  grid  spacing  remains 
constant  however,  the  x-direction  spacing  between  points  varies  with 
latitude.  The  x-direction  grid  spacing  in  Fig.  B2  was  computed  by 
taking  the  cosine  of  the  station  (p)  latitude  times  the  distance  in 
the  y-directian.  In  the  worse  case  scenario,  where  the  station 
latitude  is  the  same  as  the  latitude  of  grid  points  three  and  four  or 
grid  paints  two  and  one,  the  distance  a  would  be  approximately  1.4  km 
too  long  or  short.  Therefore,  the  approximation  of  a  square  grid  in 
the  bilinear  interpolation  formula  does  not  introduce  serious  error. 
The  computer  code  also  allows  use  of  observation  locations  outside 
the  grid  which  permits  the  use  of  exterior  grid  points  in  the 
calculations  discussed  earlier. 

Each  field  fitted  to  the  grid  was  done  using  two  scans,  the  first 
scan  used  a  radius  of  influence  of  160  km  and  the  second  scan  had  a 
radius  of  influence  of  60  km.  The  requirement  of  a  least  two 
stations  to  be  part  of  the  correction  was  imposed  on  the  first  scan 
only.  If  two  stations  were  not  within  the  radius  of  influence  during 
the  first  scan,  the  radius  of  influence  was  Increased  by  the  distance 
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between  grid  paints  in  the  y-directian  (111.  19  km)  until  the 
requirement  was  fulfilled.  The  two  scan  approach  had  a  smoothing 
effect  on  the  field  fitted  to  the  grid  point  on  the  first  pass  and 
smaller  scale  variations  within  the  field  were  introduced  in  the 
second  pass  by  the  use  of  a  much  smaller  radius  of  influence. 

However,  a  distance  weighting  scheme  is  not  the  best  objective 
analysis  procedure  according  to  Thiebaux  and  Pedder  (1987).  They 
noted  that  the  use  of  a  distance  weighting  scheme  can  lead  to  large 
biases  in  the  grid  point  data,  due  to  no  consideration  being  given  to 
the  observation  locations  relative  to  one  another.  For  example,  if 
two  stations  are  within  the  radius  of  influence  for  a  particular  grid 
point  and  they  are  both  located  in  approximately  the  same  direction 
and  distance  from  the  grid  point,  the  correction  by  each  observation 
location  is  given  the  same  weight  which  may  produce  a  somewhat 
erroneous  value  for  the  grid  point. 

APPENDIX  C  -  Computation  of  Station  Pressures,  Virtual  Temperatures, 
Gravity  Variations,  and  Data  Smoothing 

In  order  to  increase  the  number  of  observations  available  for 
each  hour,  missing  temperatures,  dewpoints,  and  altimeter  settings 
were  replaced  by  taking  the  average  value  of  the  hour  before  and 
after  the  hour  under  consideration.  If  the  respective  element  was 
missing  either  the  hour  before  or  after,  no  value  for  the  missing 
element  was  calculated.  Altimeter  settings  were  converted  to  station 
pressures  (P^y)  using 


FiS- 


Bl.  Curve  of  the  weighting  function  W  vs.  distance  d 
(After  Cressnan,  1959). 


Fig.  B2 . 


Illustration  of  the  values  used  in  the  bilinear 
interpo lat ion  formula.  (After  Inman,  1970). 
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PST  =  [  An  -  CPS  L/To>2el 


(Cl) 


where  Pe  Is  the  standard  atmospheric  sea  level  pressure  (1013.25  mb), 
A  is  the  altimeter  setting  (mb),  L e  is  the  standard  atmospheric  lapse 
rate  (6.5°  Kknf^)  ,  Te  is  the  standard  atmospheric  sea  level 
temperature  (288°K) ,  Z is  the  elevation  of  the  station  above  sea 
level  (m) ,  and  n  =  Ld  .19.  Virtual  temperature  (Tv)  was 

computed  from  the  approximation 


Tv  =  TCI  +  0. 61w),  (C2 ) 

a*  f 

where  w  is  the  actual  mixing  ratio  (Kg  Kg  ).  The  actual  mixing  ratio 
may  be  determined  by 


w  =  0.622e(PsT  -  e) 


(C3 ) 


where  e  is  the  actual  vapor  pressure  (mb).  The  actual  vapor  pressure 
was  determined  from  the  method  used  by  Koehler  (1979). 


e  *  esexp[-(L/Rv>(l/Td  -  lxT)J, 


where  e^  is  the  saturation  vapor  pressure  (mb),  L  is  the  latent  heat 
of  evaporation  (mas’*)  ,  Rv  is  the  water  vapor  gas  constant  (m*sa#K  ^  ), 
Td  is  the  dewpoint  temperature  (*K),  T  is  the  temperature,  and  a 
value  of  5430®K  was  used  for  the  ratio  L  Rv  .  The  saturated  vapor 
pressure  was  calculated  from  the  approximation  given  by  Tetens  (1930). 
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es  =6.11  X  lQCaT/b+T> 


(C5 ) 


where  the  values  used  for  the  constants  Ca  =  7.5,  b  =  237.5)  are  for 
over  a  water  surface  and  T  Is  the  temperature  (®C).  Variations  In 
gravity  (g>  due  to  elevation  and  latitude  of  the  observation  location 
was  computed  from  the  following  (Haltiner  and  Martin,  1957): 


S  =  980 . 62(1—0 . 0025 9cos20) (1-3 . 14  X  10~7z>, 


(C6  ) 


where  is  the  latitude  and  z  is  the  elevation  (m)  .  The  five-point 
formula  from  Cressman  (1959)  used  to  smooth  station  elevations  and 
the  total  f rontogenet ica 1  farcing  by  the  actual,  ageostrophic ,  and 
isal lobar ic  winds  was 


A  =  0.5A  +  1/8ZA, 


<C7) 


where  2a  is  the  sum  of  the  values  of  A  at  the  four  nearest  grid 
points. 

APPENDIX  D  -  Sensistivity  Study  on  Frontogenet ica 1  Forcing 

Winds  in  the  area  of  interest  were  measured  at  the  standard 
height  of  3,0  m  (10  feet)  above  the  surface.  Reported  wind 
directions  and  speeds  represent  a  two  minute  average  for  all 
nonmilitary  locations  (one-minute  for  Air  Force  observations).  The 
use  of  measurements  taken  as  much  apart  as  20  minutes  before  or  after 
as  being  representative  of  the  actual  surface  wind  over  the  region 


100 


opens  questions  to  the  validity  of  the  resultant  f rontogent ica 1 
forcing  for  the  actual  wind.  To  obtain  a  qualitative  measure  of 
possible  error,  (wind  direction  and  wind  speed)  50  randomly  selected 
stations  from  1900  G  M. T.  on  the  31  Jan  89  case  study  were  changed  by 
either  subtraction  or  addition  of  10°  and  2  kts,  respectively. 

After  the  results  were  obtained  for-  the  altered  conditions  given 
above,  the  same  50  observation  locations  wind  direction  and  wind 
speed  were  altered  by  random  subtraction  or  addition  of  20®and  4  kts, 
respectively.  Observation  locations  and  the  changed  wind  velocities 
are  listed  in  table  Dl.  A  count  of  the  missing  and  altered  stations 
in  table  1)1  yields  only  47  stations;  this  is  due  to  three  numbers 
being  repeated  by  the  random  number  generator  used.  Subtraction  of 
the  missing  stations  leaves  a  total  of  40  changed  wind  velocities. 
During  this  hour  293  stations  reported  wind  velocities,  hence  the  40 
changed  stations  represents  almost  14%  of  the  total  reporting 
stations.  Locations  of  the  40  stations  are  depicted  in  Fig.  Dl. 

Frontogont ica 1  forcing  by  the  actual  wind,  actual  wind  with  40 
stations  altered  by  10°  and  2  kts,  actual  wind  with  40  stations 
altered  by  20°and  4  kts  are  displayed  in  Figs.  D2-D4.  Perusal  of 
these  figtires  reveals  little  change  in  the  f  rontogenet  ical  forcing 
pattern  observed  for  the  actual  wind.  Some  magnitude  changes  are 
noted,  but  given  the  nonlinear  terms  in  Eq.  (24)  the  results  were 
better  than  anticipated.  The  smoothing  effect  of  the  Cressman  (1959) 
objective  analysis  scheme  undoubtedly  was  the  major  reason  why  more 
dramatic  changes  were  not  noted.  These  results  imply  some  error  due 
to  different  measurement  times  and  short  measurement  times  for  the 


actual  wind. 


However,  the  error  would  appear  to  be  manifested  in  the 


magnitude  and  not  in  the  horizontal  pattern  of  the  f rontogenet ica 1 
farcing. 
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Fig-  D1 .  Location  of  the  40  stations  where  the  actual  winds  were 
altered.  Stations  listed  In  Table  D1 .  are  circled. 
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Fig.  D2 .  Frontogenet lea  1  forcing  due  to  the  actual  surface  winds  at 
1900  G.  M.  T.  31  January  1989  (X  103®K  knT*hr' 1  > .  Positive 
(negative)  values  represent  f rontogenesis  (frontolysis) . 
Surface  frontal  position  is  superimposed. 


REFERENCES 


105 


Bellamy,  J.  C. ,  1945:  The  use  of  pressure  altitude  and  altimeter 

corrections  in  meteorology.  J .  Meteor .  .  2,  1-79. 

Bergeron,  T.  1937:  On  the  physics  of  fronts.  Bull.  Amer.  Meteor.  Soc . , 
18,  265-275. 

Bond,  N  A. ,  and  R.  G.  Fleagle,  1985:  Structure  of  a  cold  front  over 
the  ocean.  Quart.  J.  R.  Meteor.  Soc. ,  111,  739-759. 

Bonner,  W.  D. ,  and  J.  Paegle,  1970:  Diurnal  variations  in  boundary 

layer  winds  over  the  south-central  United  States  in  summer.  Mon . 
Vea .  Rev. ,  98,  735-744. 

Browning,  K.  A. ,  and  G.  A.  Monk,  1982:  A  simple  model  for  the  synoptic 
analysis  of  cold  fronts.  Quart.  J,  R.  Meteor.  Soc.  ,  108,  435-452. 

Brundidge,  K.  C. ,  1965:  The  wind  and  temperature  structure  of 

nocturnal  cold  fronts  in  the  first  1420  feet.  Mon.  Vea.  Rev.  , 

93,  587-603. 

Brunt,  D. ,  and  C.  K.  M.  Douglas,  1928:  The  modification  of  the 

strophic  balance  for  changing  pressure  distribution,  and  its 
effect  on  rainfall.  Mem.  Roy.  Meteor.  Soc. ,  3,  29-51. 

Cressman,  G.  P.,  1959:  An  operational  objective  analysis  system.  Mon. 

Vea.  Rev. ,  87,  367-374. 

Ualtiner,  G.  J. ,  and  F.  L.  Martin,  1957:  Dynamical  and  Physical 
Meteorology.  New  York,  McGraw-Hill,  470pp. 

Hobbs,  P.  V. ,  and  P.  0.  V.  Pearson,  1982:  The  mesoscale  and  microscale 
structure  and  organization  of  clouds  and  precipitation  in 
midlatitude  cyclones.  Part  V:  The  struture  of  narrow  cold-frontal 
rainbands.  J .  Atmos .  Sc  1 .  ,  39,  280-295. 

Hoskins,  B.  J. ,  and  F.  P.  Bretherton,  1972:  Atmospheric  f rontogenesis 
models:  Mathematical  formulation  and  solution.  .J  .  Atmos ■  Sc  1  ■  ,  29, 

11-37, 

Inman,  R.  L. ,  1970:  Papers  on  operational  objective  analysis  schemes 

at  the  National  Severe  Storms  Forecast  Center.  Tech.  Memo 
ERTI.M  NSSL  51.  ,  91  pp. 

Jascourt,  S.  D. ,  S.  S.  Lindstrora,  C.  J.  Seman,  and  D.  D.  Houghton, 

1988:  An  observation  of  banded  convective  development  in  the 
presence  of  weak  symmetric  stability.  Mon ,  Vea .  Rev ■ ,  116, 

175-191. 

Koehler,  T.  L. ,  1979:  A  case  study  of  height  and  temperature  analysis 

derived  from  Nimbus-6  satellite  soundings  on  a  fine  mesh  model 
grid.  Ph.D.  Dissertation,  Department  of  Meteorology,  University 
of  Vi sconsin-Mad Ison ,  185  pp. 


106 


McDonell,  J.  E. ,  1967:  A  summary  of  the  first-guess  fields  used  at  the 

National  Meteorological  Center.  Tech.  Memo  WBTMNM3-38. ,  U.  S. 
Weather  Bureau,  17  pp. 

Moore,  J.  T. ,  and  P.  D.  Blakley,  1988:  The  role  of  f rontogenet ical 
forcing  and  conditional  symmetric  instability  in  the  midwest 
snowstorm  of  30-31  January  1982.  Mon.  We a.  Rev. ,  116,  2155-2171. 

Petterson,  S. ,  1956:  Weather  Analysis  and  Forecasting,  Vol  1. 

McGraw-Hill,  428  pp. 

Reeder,  M.  J . ,  1986:  The  interaction  of  a  cold  front  with  a  prefrontal 

thermodynamically  well-mixed  boundary  layer.  Aust.  Meteor.  Mag. . 
34,  137-148. 

- ,  and  R.  K.  Smith,  1986:  A  comparison  between  f rontogenes is  in 

the  two-dimensional  Eady  model  of  baroclonic  instability  and 
summertime  cold  fronts  in  the  Australian  region.  Quart .  J .  Rov . 
Meteor.  Soc ■ ,  112,  293-313. 

- ,  and  - ,  1987:  A  study  of  frontal  dynamics  with  application 

to  the  Australian  summertime  "cool  change.”  J ■  Atmos.  Scl . ,  44, 

687-705. 

Sangster,  W. E. ,  1960:  A  method  of  representing  the  horizontal  pressure 

gradient  force  without  reduction  of  station  pressures  to  sea 
level.  J .  Meteor. ,  17,  166-176. 

- ,  1967:  Diurnal  surface  geostrophic  wind  variations  over  the 

Great  Plains.  Preprints,  Fifth  Conf.  on  Severe  Local  Storms,  St. 
Louis  ,  MO,  146-153.  Amer.  Meteor.  Soc. 

- ,  1987:  An  improved  technique  for  computing  the  horizontal 

pressure-gradient  force  at  the  Earth's  surface.  Mon.  Wea .  Rev., 

115,  1358-1369. 

Saucier,  W.  J.  ,  1955:  Principles  of  Meteorological  Analysis.  Chicago 

Press,  438  pp. 

Smith,  R.  K.  ,  and  M.  J.  Reeder,  1988:  On  the  movement  and  low  level 
structure  of  cold  fronts.  Mon .  Wea .  Rev ■ ,  116,  1927-1944. 

Tetans,  O.  1930:  Z.  Geophysics. ,  6. 

Thiebaux,  H.  J.,  and  M.  A.  Pedder,  1987:  Spatial  Objective  Analysis: 
with  applications  In  atmospheric  science.  Academic  Press,  299  pp. 

Young,  J.  A. ,  1973:  A  theory  for  isallobaric  air  flow  in  the  planetary 

boundary  layer.  J .  Atmos.  Sci , ,  30,  1584-1592. 


